SPRINGER BRIEFS IN EARTH SYSTEM SCIENCES 


Marine Carbon 
Biogeochemistry 
A Primer for Earth 
= System Scientists 


Z) Springer Open 


SpringerBriefs in Earth System Sciences 


Series editors 


Gerrit Lohmann, Universität Bremen, Bremen, Germany 

Lawrence A. Mysak, Department of Atmospheric and Oceanic Sciences, McGill 
University, Montreal, QC, Canada 

Justus Notholt, Institute of Environmental Physics, University of Bremen, Bremen, 
Germany 

Jorge Rabassa, Laboratorio de Geomorfologia y Cuaternario, CADIC-CONICET, 
Ushuaia, Tierra del Fuego, Argentina 

Vikram Unnithan, Department of Earth and Space Sciences, Jacobs University 
Bremen, Bremen, Germany 


SpringerBriefs in Earth System Sciences present concise summaries of cutting-edge 
research and practical applications. The series focuses on interdisciplinary research 
linking the lithosphere, atmosphere, biosphere, cryosphere, and hydrosphere 
building the system earth. It publishes peer-reviewed monographs under the 
editorial supervision of an international advisory board with the aim to publish 8 to 
12 weeks after acceptance. Featuring compact volumes of 50 to 125 pages (approx. 
20,000—70,000 words), the series covers a range of content from professional to 
academic such as: 


A timely reports of state-of-the art analytical techniques 
bridges between new research results 

snapshots of hot and/or emerging topics 

literature reviews 

in-depth case studies 


Briefs are published as part of Springer’s eBook collection, with millions of users 
worldwide. In addition, Briefs are available for individual print and electronic 
purchase. 

Briefs are characterized by fast, global electronic dissemination, standard 
publishing contracts, easy-to-use manuscript preparation and formatting guidelines, 
and expedited production schedules. 

Both solicited and unsolicited manuscripts are considered for publication in this 
series. 


More information about this series at http://www.springer.com/series/10032 


Jack J. Middelburg 


Marine Carbon 
Biogeochemistry 


A Primer for Earth System Scientists 


A Springer Open 


Jack J. Middelburg 
Department of Earth Sciences 
Utrecht University 

Utrecht, The Netherlands 


ISSN 2191-589X ISSN 2191-5903 (electronic) 
SpringerBriefs in Earth System Sciences 
ISBN 978-3-030-10821-2 ISBN 978-3-030-10822-9 (eBook) 


https://doi.org/10.1007/978-3-030-10822-9 
Library of Congress Control Number: 2018965889 


© The Editor(s) (if applicable) and The Author(s) 2019. This book is an open access publication. 
Open Access This book is licensed under the terms of the Creative Commons Attribution 4.0 
International License (http://creativecommons.org/licenses/by/4.0/), which permits use, sharing, adap- 
tation, distribution and reproduction in any medium or format, as long as you give appropriate credit to 
the original author(s) and the source, provide a link to the Creative Commons licence and indicate if 
changes were made. 

The images or other third party material in this book are included in the book’s Creative Commons 
licence, unless indicated otherwise in a credit line to the material. If material is not included in the book’s 
Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the 
permitted use, you will need to obtain permission directly from the copyright holder. 

The use of general descriptive names, registered names, trademarks, service marks, etc. in this publi- 
cation does not imply, even in the absence of a specific statement, that such names are exempt from the 
relevant protective laws and regulations and therefore free for general use. 

The publisher, the authors, and the editors are safe to assume that the advice and information in this 
book are believed to be true and accurate at the date of publication. Neither the publisher nor the 
authors or the editors give a warranty, express or implied, with respect to the material contained herein or 
for any errors or omissions that may have been made. The publisher remains neutral with regard to 
jurisdictional claims in published maps and institutional affiliations. 


This Springer imprint is published by the registered company Springer Nature Switzerland AG 
The registered company address is: Gewerbestrasse 11, 6330 Cham, Switzerland 


Preface 


Biogeochemistry, a branch of Earth System Sciences, focusses on the two-way 
interactions between organisms and their environment, including the cycling of 
energy and elements and the functioning of organisms and ecosystems. To this end, 
physical, chemical, biological and geological processes are studied using field 
observations, experiments, modelling and theory. The discipline of biogeochem- 
istry has grown to such an extent that sub-disciplines have emerged. Consequently, 
producing a single comprehensive textbook covering all aspects, e.g., terrestrial, 
freshwater and marine domains, biogeochemical cycles and budgets of the major 
biological relevant elements, reconstruction of biogeochemical cycles in the past, 
earth system modelling, microbiological, organic and inorganic geochemical 
methods, theory and models, has become unworkable. 

This book provides a concise treatment of the main concepts in ocean carbon 
cycling research. It focusses on marine biogeochemical processes impacting the 
cycling of particulate carbon, in particular organic carbon. Other biogeochemical 
processes impacting nitrogen, phosphorus, sulphur, etc., and the identity of the 
organisms involved are only covered where needed to understand carbon biogeo- 
chemistry. Moreover, chemical and biological processes relevant to carbon cycling 
are central, i.e. for physical processes, the reader might consult the excellent ocean 
biogeochemical dynamics textbooks of Sarmiento and Gruber (2006; Princeton 
University Press) and Williams and Follows (2011; Cambridge University Press). 
My text aims to provide graduate students in marine and earth sciences a conceptual 
understanding of ocean carbon biogeochemistry, so that they are better equipped to 
read palaeorecords, can improve carbon biogeochemical models and generate more 
accurate projections of the functioning of the future ocean. Because the book is 
targeted at students having a background in environmental and earth sciences, some 
basic biological concepts are explained. Some basic understanding of calculus is 
expected. Simple mathematical models are used to highlight the most important 
factors governing carbon cycling in the ocean. The material here is based on a 
selection of lectures in my Utrecht University master course on Microbes and 
Biogeochemical Cycles. 

This first draft of this book was written during a three-month sabbatical stay at 
Department of Geosciences, Princeton University (April-June 2018). I thank Bess 
Ward, chair of that department, for providing a desk and a stimulating environment. 
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This sabbatical stay was supported by a travel grant from the Netherlands Earth 
System Science Centre. I thank Bernie Boudreau for carefully scrutinizing the 
initial draft, Mathilde Hagens and Karline Soetaert for feedback on Chap. 5 and 
Anna de Kluijver for remarks on Chap. 6. Ton Markus improved my draft figures. 
Finally, I thank my wife and publisher Petra van Steenbergen. 


Utrecht, The Netherlands Jack J. Middelburg 
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Symbols 


Biomass of phytoplankton or buffer value (Chap. 5) 

Diffusion coefficient (area time |); with D,: diffusion of solutes in 
sediments, D,: particle mixing in sediments 

Radiant energy (mol quanta area ' time ‘) 

Activation energy (J mol ') 

Flux of material (mol/gr area ' time ') 

Quantity of organic carbon (mol/gr C per gr sediment, area or volume) 
First-order rate/decay constant (time ') 

Light extinction coefficient (length ') 

Half-saturation constant in Monod-type equation; Kg: light saturation 
parameter; K,,: growth (nutrient) half-saturation constant 

Equilibrium constants (x = w, H, 1, 2) that depend on temperature, pressure 
and solution composition 

Eddy-diffusion (mixing) coefficient in water column (area time ') 
Production (mol/gr volume ' time‘) 

Increase in rate for 10 °C increase in T 

Cellular quota in Droop equation 

First-order rate constant for phytoplankton (time™') 

Zero-order production or consumption term (mol/gr volume ' time ') 
Time 

Temperature (°C or K) 

Particle settling in water column or sediment accumulation rate (length 
time ') 

Depth in sediment (length) 

Depth in water column (length) 

Euphotic zone depth (length) 

Compensation depth (length) where phytoplankton growth and respiration 
are equal 

Critical depth (length) where phytoplankton production balances losses 
Buffer value in terms of proton concentration 

Solute transfer coefficient at seafloor 


x< 


<D DES 


Porosity in sediment 

Maximum growth rate (time!) 

Maximum nutrient uptake 

Dry density of particle (gr volume ') 

Carbon dioxide generated per unit carbonate precipitated 


Symbols 


Introduction 


The name biogeochemistry implies that it is a discipline integrating data, knowl- 
edge, concepts and theory from biology, geosciences and chemistry. Biogeo- 
chemists extensively use approaches from a wide range of disciplines, including 
physical, chemical and biological oceanography, limnology, atmospheric sciences, 
ecology and microbiology, civil and environmental engineering, soil science and 
geochemistry. This diversity in scientific backgrounds stimulates cross-fertilization 
and research creativity, which are needed to elucidate the reciprocal relationships 
between living organisms and their environment at multiple scales during times of 
global change. Biogeochemistry aims to provide a holistic picture of natural 
ecosystem functioning. The challenge is to identify the right level of detail needed 
to understand the dynamics of elemental cycles and the functioning of biological 
communities. This implies that single-cell organism level studies and molecular 
orbital calculations of chemical reactions require upscaling to the appropriate 
temporal and spatial scale (often involving first-principle physics based models) to 
understand how natural ecosystems deal with perturbations and how life has shaped 
our planet. 

Although biogeochemistry developed as a full discipline in the mid-1980s with 
the launch of the international geosphere-biosphere program (IGBP, 1987) and the 
journals Biogeochemistry (1984) and Global Biogeochemical Cycles (1987), its 
roots can be traced back to early scientists documenting how living organisms 
transformed chemical substances, such as oxygen production during photosynthesis 
(Priestly, 1733—1804), phosphorus in organisms’ tissues (Lavoisier, 1743-1789) 
and nitrogen fixation by bacteria (Beijerinck, 1851-1931). Naturalist and 
avant-la-lettre multidisciplinary scientists, such as Alexander von Humboldt (1769- 
1859). Charles Darwin (1808-1882) and Alfred Lotka (1880-1949), pioneered 
what we would recognize as biogeosciences in the 21st century. Darwins’ studies of 
atmospheric deposition, bioturbation and formation and sustenance of coral reefs 
are still key areas in modern biogeochemistry. The tight relationship between living 
organisms and their environment figured prominently in Lotka’s book “Elements of 
Physical Biology” (1925): “Jt is not so much the organism or the species that 
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evolves, but the entire system, species and environment. The two are inseparable.” 
This concept that organisms shape the environment and govern elemental cycles on 
Earth underlies the biosphere concept of Vladimir Vernadsky (1863-1945), a 
geochemist and mineralogist, often considered the founder of biogeochemistry. 
G. Evelyn Hutchinson (1903-1991) was instrumental in establishing biogeo- 
chemical, whole-system approaches to study lakes. Alfred Redfield (1890-1983) 
discovered that nitrogen to phosphorus ratios of phytoplankton in seawater are 
constant and similar to dissolved ratios, implying co-evolution of the environment 
and organisms living in it. His seminal 1958 article started as follows “It is a 
recognized principle of ecology that the interaction of organisms and environment 
are reciprocal. The environment not only determines the conditions under which 
life exists, but the organisms influence the conditions prevailing in their environ- 
ment” (Redfield 1958). The latter was articulated in the Gaia hypothesis of Love- 
lock (1972): The Earth became and is maintained habitable because of multiple 
feedback mechanisms involving organisms. For instance, biologically mediated 
weathering of rocks removes carbon dioxide from the atmosphere and generates 
bicarbonate and cations that eventually arrive in the ocean, where calcifiers produce 
the minerals calcite and aragonite and release carbon dioxide back to the 
atmosphere. 

The above one-paragraph summary of the history of biogeochemistry does not 
mean that it was a linear or smooth process. While the early pioneers (before the 
second world war) were not hindered much by disciplinary boundaries between 
physics, biology, chemistry and earth sciences, the exponential growth of scientific 
knowledge and the consequent specialization and success of reductionism to 
advance science, had led to an under appreciation of holistic approaches crossing 
disciplinary boundaries during the period 1945-1990. Addressing holistic research 
questions may require development of new concepts and methods, but often 
involves application and combination of well-established theory or methods from 
multiple disciplines. The latter implies finding the optimal balance between biology, 
chemistry and physics to advance our understanding of biogeochemical processes. 
For instance, all biogeochemical models have to trade-off spatial resolution in the 
physical domain with the number of chemical elements/compounds and the 
diversity of organisms to be included. Ignoring spatial dimensions and hetero- 
geneity through the use of box models may seem highly simplistic to a physical 
oceanographer, but may be sufficient to obtain first-order understanding of ele- 
mental cycling. Similarly, organic carbon flows can be investigated via study of the 
organisms involved, the composition of the organic matter or by quantifying the 
rates of transformation, without considering the identity of the organisms involved. 
Each disciplinary approach has its strengths and weaknesses, and they are unfor- 
tunately not always internally consistent. However, this confrontation of different 
disciplinary concepts has advanced our understanding (Middelburg 2018). In the 
next section, we will discuss why many geochemists embraced biogeochemistry. 
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1.1 From Geochemistry and Microbial Ecology 
to Biogeochemistry 


Geochemistry is a branch of earth sciences that applies chemical tools and theory to 
study earth materials (minerals, rocks, sediments and water) to advance under- 
standing of the Earth and its components. While early studies focused on the 
distribution of elements and minerals using tools from analytical chemistry, the next 
step involved the use of chemical thermodynamics to explain and predict the 
occurrence and assemblages of minerals in sediments and rocks. The thermody- 
namic approach was and is very powerful in high-temperature systems (igneous 
rocks, volcanism, metamorphism, hydrothermal vents), but it was less successful in 
predicting geochemical processes at the earth surface. Geochemists studying earth 
surface processes soon realized that predictions based on thermodynamics, i.e. the 
Gibbs free energy change of a reaction, provided a necessary condition whether a 
certain reaction could take place, but not a sufficient constraint whether it would 
take place because of kinetics and biology. 

Realizing the limitations of the thermodynamic approach, the field of geo- 
chemical kinetics developed from the 1980s onwards (Lasaga 1998). Much pro- 
gress was made studying mineral precipitation and dissolution kinetics as a function 
of solution composition (e.g. pH) and environmental conditions (e.g., temperature). 
These laboratory studies were done under well constrained conditions and in the 
absence of living organisms. However, application of these experimentally deter- 
mined kinetic parameters to natural systems revealed that chemical kinetics often 
could not explain the differences between predictions based on chemical thermo- 
dynamics and kinetics, and observations in natural systems. These unfortunate 
discrepancies were attributed to the black box ‘biology’ or ‘bugs’. 

Before the molecular biology revolution, microbial ecology was severely 
method limited. Samples from the field were investigated using microscopy and 
total counts of bacteria were reported. Microbiologists were isolating a biased 
subset of microbes from their environment and studying their metabolic capabilities 
in the laboratory. To investigate whether these microbial processes occur in nature, 
microbial ecologists developed isotope and micro-sensor techniques to quantify 
rates of metabolism in natural environments (e.g., oxygen production or con- 
sumption, carbon fixation, sulfate reduction). These microbial transformation rates 
were of interest to geochemists because they represented the actual reaction rates, 
rather than the ones predicted from geochemical kinetics. Microbial ecologists and 
geochemists started to collaborate systematically and a new discipline emerged in 
which cross-fertilization of concepts, approaches and methods stimulated not only 
research questions at the interface but also in the respective disciplines. Stable 
isotope and organic geochemical biomarker techniques and detailed knowledge on 
mineral phases have enriched geomicrobiology, while knowledge on microbes and 
their capabilities and activities has advanced the understanding of elemental 
cycling. This integration of microbial ecology and geochemistry has evolved well 
regarding tools (e.g., the use of compound-specific isotope analysis and nanoSIMS 
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in microbial ecology for identity-activity measurements), but less so in terms of 
concepts and theoretical development. Moreover, there is more to biology than 
microbiology. Animals and plants have a major impact on biogeochemical cycles, 
not only via their metabolic activities (primary production, nutrient uptake, respi- 
ration), but also via their direct impact on microbes (grazing, predation) and their 
indirect impact via the environment (ecosystem engineering: e.g., bioturbation, soil 
formation). This additional macrobiological component of biogeochemistry is 
increasingly being recognized (Middelburg 2018). 


1.2 Focus on Carbon Processing in the Sea 


This book focuses on biogeochemical processes relevant to carbon and aims to 
provide the reader (graduate students and researchers) with insight into the func- 
tioning of marine ecosystems. A carbon centric approach has been adopted, but other 
elements are included where relevant or needed; the biogeochemical cycles of 
nitrogen, phosphorus, iron and sulfur are not discussed in detail. Furthermore, the 
organisms involved in carbon cycling are not discussed in detail for two reasons. 
First, this book focuses on concepts and the exact identity of the organisms involved 
or the systems (open ocean, coastal, lake) is then less relevant. Secondly, our 
knowledge of the link between organism identity and activity in natural environments 
is limited. For instance, primary production rates are often quantified and phyto- 
plankton community composition is characterized as well, but their relationship is 
poorly known. The extent of particle mixing by animals in sediments can be quan- 
tified and the benthic community composition can be described, but the contribution 
of individual species to particle mixing cannot be estimated in a simple manner. 
The following chapters will respectively deal with production (Chap. 2) and 
consumption (Chap. 3) of organic carbon in the water column, the processing of 
organic carbon at the seafloor (Chap. 4), the impact of biogeochemical processes on 
inorganic carbon dynamics (Chap. 5), and the composition of organic matter 
(Chap. 6). The carbon cycle is covered using concepts, approaches and theories 
from different subdisciplines within ecology (phycologists, microbial ecologists and 
benthic ecologists) and geochemistry (inorganic and organics) and crosses the 
divides between pelagic and benthic systems, and coastal and open ocean. The book 
aims to provide the reader with enhanced insight via the use of very simple, generic 
mathematical models, such as the one presented in Box 1.1. Because of our focus 
on concepts, in particular the biological processes involved, there will be little 
attention to biogeochemical budgets and the role of large-scale physical processes 
in the ocean (Sarmiento and Gruber 2006; Williams and Follows 2011). Accurate 
carbon budgets are essential for a first-order understanding of biogeochemical 
cycles, but it is important to understand the mechanisms involved before adequate 
projections can be made for the functioning of System Earth and its ecosystems in 
times of change. To set the stage for a detailed presentation of biogeochemical 
processes, we first introduce a simple organic carbon budget for the ocean. 
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1.3 A 101 Budget for Organic Carbon in the Ocean 


Establishing carbon budgets in the ocean, in particular during the Anthropocene, is 
a far from trivial task, involving assimilation of synoptic remote sensing and sparse 
and scarce field observations with deep insight and numerical modelling of the 
transport and reaction processes in the ocean. The important processes and thus 
flows of carbon in the ocean are related to primary production, export of organic 
carbon from the surface layer to ocean interior, deposition of organic carbon at the 
seafloor and organic carbon burial in sediments. Accepting 25% uncertainty, these 
numbers are well constrained at 50 Pg C y' (1 Pg or 1 Gt is 10° gr) for net 
primary production, 10 Pg C y | for export production, 2 Pg C y~! for carbon 
deposition at the seafloor and 0.2 Pg C y ' for organic carbon burial (Fig. 1.1). 
Although no detailed, closed complete carbon budgets will be presented, estimates 
for individual processes, including gross primary production, chemoautotrophy and 
coastal processes, are presented in the following chapters. However, the 
50-10-2-0.2 rule for carbon produced, transferred to the ocean interior, deposited at 


Fig. 1.1 Simplified budget 
of carbon flows in the ocean. 
Each year net phytoplankton 


production is about 50 Pg C Euphotic -1 

(1 Pg = 1 Gt = 10" g), 10 Pg zone 207100. m 
is exported to the ocean 

interior, the other 40 Pg is 50 Produced 


respired in the euphotic zone. 
Organic carbon degradation 
continues while particles 
settle through the ocean 10 Exported 
interior and only 2 Pg 
eventually arrives at the 
seafloor, the other 8 Pg is 
respired in the dark ocean. In 
sediments, the time scale 
available for degradation 
increases order of magnitude 
with the result that 90% of the 
organic carbon delivered is 
degraded and only 0.2 Pg C 
yr | is eventually buried and 
transferred from the biosphere 
to the geosphere 


Ocean 


Interior 3-5 km 


2 Deposited 


Sediments 10-30 cm 
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the seafloor and preserved in sediments, respectively, can easily be remembered and 
should be kept in mind when reading the details of carbon processing in the 
remaining of this book. 


Box 1.1: A simple mathematical model for reaction and transport 

In multiple chapters, we will make use of a very simple mathematical model 
in which the change in C (concentration, biomass) is due to the balance 
between diffusion (eddy Kz, molecular D), advection (sediment accretion 
particle/phytoplankton settling, w) and net effects of reactions (production 
and consumption). The basic equation is: 


oc OC oC 
ot a Ge WOK Be 


aC 2C is 
at Ox? 
the spatial change in transport due to diffusion with diffusion coefficient D 
(m? s"'), w as is the spatial change in transport due to water flow or particle 
settling with velocity w (m s '), positive downwards, —kC is the con- 
sumption of substance C via a first order reaction with reactivity constant k 
(s7!) and Ro is a zero-order production term (that is, the substance C has no 
impact on the magnitude of this rate). 

This equation is based on spatially uniform mixing and settling rates and 
reactivity (ie. D, w and k are constant). Moreover, we consider only 


steady-state conditions, i.e. there is no dependence on time. This simplifies 


where is the change in concentration (mol m °) with time (t, s), D 


the math: the partial differential equation (%) becomes an ordinary differ- 
a s dC). 
ential equation (aC): 
D ae ae kC+Ro = 0 
w t = 
dx? dx : 


If we first consider the situation without zero-order production or consump- 
tion (i.e. Ro = 0), the general solution is: 


C = Ae™ + Be™ 


w — y w? + 4kD w + V w2 + 4kD 
where « = — p p= — 


and A and B are integration constant depending on the boundary conditions. 


The number of integration constants sets the number of boundary conditions 


required. We will use models for the semi-infinite domain: i.e., if x — œ then 


the gradient in C disappears (JE = 0). Since all terms in B are positive, the 
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second term becomes infinite and the integration constant B must thus be zero 
for this boundary condition. 

For the upper boundary condition, we will explore two types: a fixed 
concentration and a fixed flux condition. If we know C = Co at depth x = 0, 
then A is Co and the solution is: 


(C= Coe™. 


Sometimes we know the external flux (F) of C, then we have to balance the 
flux at the interface at x = 0, e.g.: 


dC 
F=-D=—| wees 
dx x=0 


Next, we take the derivative of the remaining first-term of the general solution 
(Ae™), to arrive at: 


F = —DaAe” + wAe” 
Since e? = 1, A= aw and the solution is: 


F 
C = — e”. 
Duw 


In some systems, transport is dominated by diffusion (e.g. molecular diffusion 


of oxygen in pore water, eddy diffusion of solutes and particles in water) and 
the advection term (w) can be ignored. The basic solutions given above 


remain but now & = — iS and the pre-exponential term for the constant flux 


upper boundary becomes — ie In other systems transport is dominated by the 


advection term (e.g. settling particles in the water column) and then « = — k 


Ww 
and the flux upper boundary condition becomes BI 

The above solutions are valid in the case that only first-order reaction 
occurs. The presence of zero-order reactions results in different solutions and 
these will be presented in the text where relevant. Similarly, the solutions 
presented are only valid if D, w and k are uniform with depth. In Chap. 3 we 
present an advection-first order degradation model in which we vary w and k 
with depth. Although user-friendly packages and accessible textbooks are 
available for numerical solving these and more complex equations (Boudreau 
1997; Soetaert and Herman 2009), we restrict ourselves to analytical solutions 
because the relations among D, w and k in the various applications reveal 
important insights in the various process and governing factors, and the 
reader can implement the analytical solutions for further study. 
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Primary Production: From Inorganic 2 
to Organic Carbon 


Primary production involves the formation of organic matter from inorganic carbon 
and nutrients. This requires external energy to provide the four electrons needed to 
reduce the carbon valence from four plus in inorganic carbon to near zero valence in 
organic matter. This energy can come from light or the oxidation of reduced 
compounds, and we use the terms photoautotrophy and chemo(litho)autotrophy, 
respectively. Total terrestrial and oceanic net primary production are each ~ 50-55 
Pg yr | (1 Pg = 1 Gt = 10" g; Field et al. 1998). Within the ocean, carbon fixation 
by oceanic phytoplankton (~47 Pg yr ') dominates over that by coastal phyto- 
plankton (~6.5 Pg yr '; Dunne et al. 2007), benthic algae (~ 0.32 Pg yr '; Gattuso 
et al. 2006), marine macrophytes (~1 Pg yr; Smith 1981) and chemo(litho) 
autotrophs (~0.4 and ~0.37 Pg yr ' in the water column and sediments, 
respectively; Middelburg 2011). Much of the chemolithoautrophy is based on 
energy from organic matter recycling. Since, photosynthesis by far dominates 
inorganic to organic carbon transfers, we will restrict this chapter to light driven 
primary production. 

Gross primary production refers to total carbon fixation/oxygen production, 
while net production refers to growth of primary producers and is lessened by 
respiration of the primary producer. Net primary production is available for growth 
and metabolic costs of heterotrophs, and it is the process most relevant for bio- 
geochemists and chemical oceanographers. For the time being, we present primary 
production as the formation of carbohydrates (CH2O) and ignore any complexities 
related to the formation of proteins, membranes and other cellular components 
(Chap. 6), because these require additional elements (nutrients). The overall pho- 
tosynthetic reaction is: 


CO2 + H20 + light — CH20 + O2 


© The Author(s) 2019 9 
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It starts with the absorption of light energy by photosystem II (PSII): 
2H0 + light  4H*+ +4e7 + Op 


This reaction yields energy to generate adenosine triphosphate (ATP). The oxygen 
produced originates from the water and can be considered a waste product of 
photosynthesis. The protons and electrons generated subsequently react with 
nicotinamide adenine dinucleotide phosphate (NADP*) at photosystem I (PSI): 


NADP* +H* +2e7 NADPH. 


The energies of NADPH and ATP are then used to fix and reduce CO, to form 
carbohydrate. 


CO, +4H* + 407 25C? CH O + HO 


This reaction is normally mediated by the enzyme ribulose bis-phosphate car- 
boxylase (RuBisCO). 

Primary production is at the base of all life on earth; it is thus important to 
quantify it and to understand the governing factors. We will first present, at a very 
basic level, the primary producers. This will be followed by the introduction of the 
master equation of primary production, based on laboratory studies, and then a 
discussion of its application to natural systems. 


2.1 Primary Producers 


Primary producers in the ocean vary from pm-sized phytoplankton to m-sized 
mangrove trees. Phytoplankton refers to photoautotrophs in the water that are 
transported with the currents (although they may be slowly settling). Biological 
oceanographers usually divide plankton (all organisms in the water that go with the 
current) into size classes (Table 2.1). Most phytoplankton are in the pico, nano and 
microplankton range (0.2—200 um). The prefixes pico and nano have little to do 
with their usual meaning in physics and chemistry. Their small size gives them a 
high-surface-area-to-volume ratio which is highly favourable for taking up nutrients 
from a dilute solution. Within these phytoplankton size classes there is high 
diversity in terms of species composition and ecological functioning. Both small 
cyanobacteria (Synechococcus and Procholoroccus) and very small eukaryotes 
(e.g., Chlorophytes) contribute to the picoplankton. Microflagellates from various 
phytoplankton groups (Chlorophytes, Cryptophytes, Diatoms, Haptophytes) dom- 
inate the nanoplankton and differ in many aspects (cell wall, nutrient stoichiometry, 
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Table 2.1 Plankton size classes in the ocean 


Size lass Name (example) 

<0.2 um Femtoplankton (virus) 

0.2-2 um Picoplankton (bacteria, very small eukaryotes) 
2-20 um Nanoplankton (diatoms, dinoflaggelates, protozoa) 
20-200 um Microplankton (diatoms, dinoflaggelates, protozoa) 
0.2-20 mm Mesoplankton (zooplankton) 

2-20 cm Macroplankton 


pigments, number of flagellae, life history, presence/absence of frustule). While 
phytoplankton communities can be described in terms of species, size classes or 
molecular biology data based partitioning units, they can also be divided into 
different functional types (diatoms because of Si skeleton, coccoliths with CaCO3 
skeleton, N>-fixers, etc.). Unfortunately, taxonomic, functional and size partition- 
ings among phytoplankton groups are not necessarily consistent. 

A substantial fraction of the ocean floor in the coastal domain receives enough 
light energy to sustain growth of photoautotrophs. This includes not only intertidal 
areas, but also the subtidal. Small-sized photoautotrophs (microphytobenthos, 
including diatoms and cyanobacteria) are again the dominant primary producers, 
but macroalgae, seagrass, saltmarsh plants and mangrove trees contribute as well. 
Seagrasses, saltmarsh macrophytes and mangrove trees have structural components 
and specialised organs (roots and rhizomes) to tap into nutrient resources within the 
sediments. 


2.2 The Basics (For Individuals and Populations) 


Carbon fixation by (and growth of) primary producers will be discussed based on 
the master equation of Soetaert and Herman (2009): 


P = u- B - flim (resources, conditions) (2.1) 


This master equation simply states that production (P, mol/g per unit volume per 
unit time) is proportional to the biomass (B, mol/g per unit volume) of the primary 
producer, the actor, which has an intrinsic maximum growth rate of u (time ') and 
is limited (0 < fiim < 1) by either physical conditions (e.g., temperature, turbulence) 
or resources such as light, nutrients and dissolved inorganic carbon. This equation is 
simple and generic, and we will show below how it relates to phytoplankton global 
primary production estimates using remote sensing, to expressions used in 
numerical biogeochemical models and to exponential growth in the laboratory. 
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2.2.1 Maximum Growth Rate (p) 


Consider a primary producer in an experiment supplied with all the resources it 
needs and under ideal conditions, in other words the limitation function fiim is equal 
to one and optimal growth occurs. Equation 2.1 then reduces to the change in B 
with time, or production P, is equal to p-B: 


dB 


P=—=uB 2.2 
g7" (2.2) 


This is the well-known equation for exponential growth: 


L Roel ie es 
B = Boe", or alternatively : u : ing (2.3) 
0 


where B is the biomass at times t and Bo is the initial biomass. Plotting the 
logarithm of biomass development as function of time yields then a slope corre- 
sponding to u. Sometimes data are reported as the number of cell divisions (or 
doublings) per day: u4 = } logo 5- 

Maximum growths for phytoplankton typically varies from 0.1 to 4 d', 
implying doubling times (2) of a fraction of a day to one week. Figure 2.la 


shows a typical example of exponential growth for maximum growth rates of 0.1 to 
2 d-'. Exponential growth leads to rapid depletion of substrates and after some 
time, resources become limiting and phytoplankton enters into a stationary phase 
(Fig. 2.1b). Maximum growth size depends on phytoplankton group and size 
(Fig. 2.2; Box 2.1). 
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Fig. 2.1 a The increase in biomass during exponential growth with growth rates of 0.1, 0.5, 1 and 
2 d™!. b Cell growth of the diatom Thalassiosira pseudonana is exponential (growth rate of 1.4 
d-') till nitrate is depleted and then stationary growth occurs (Data from Davidson et al. 1999) 
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Fig. 2.2 a The relationship between growth rate and phytoplankton cell size and model 
prediction (grey curve) (Ward et al. 2017). b The relationship between growth rate and 
phytoplankton cell size (solid black line). Maximum nutrient uptake and requirement per cell scale 
positively with cell size (dashed blue line), while theoretical maximum growth rates scale 
negatively (solid blue line) 


2.2.2 Temperature Effect on Primary Production 


The temperature of a system provides a strong control on the functioning of 
organisms. Growth responses of populations to temperature are usually expressed 
by thermal tolerance curves, also known as reaction norms. Starting at low tem- 
peratures, growth initially increases linearly or exponentially up to a maximum Topt 
and then typically declines relatively more rapidly: i.e. the response curve is often 
skewed to the left. In other words, phytoplankton growing near its optimum tem- 
perature is more sensitive to warming than to cooling (Fig. 2.3). 

Although populations show distinct unimodal responses to temperature, mixed 
communities, and thus ecosystems, usually exhibit a smooth, monotonical increase 
best described by an exponential (u = aet, Fig. 2.4). The thermal response can 
then be described by 


Pat 
bT opt 
= 1 — | ——— 2.4 
ae Gena oa) 


where a and b are empirical parameters describing the maximum envelope for the 
mixed community and T,,, and width describe the maximum growth rate and 
temperature range of individual populations. Eppley’s (1972) seminal work on 
temperature and phytoplankton growth in the sea reported values of 0.59 for a and 
0.0633 for b. Note that this community response provides an upper limit for 


individual species and that high growth rates for individual species trade off with 
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Fig. 2.3 Temperature response of the specific growth rates of three eukaryotic phytoplankton 
classes (Diatoms, Dinoflaggelates, Chlorophytes) and of Cyanobacteria from temperate freshwater 
and brackish environments (modified from Paerl et al. 2011). A. form = Asterionella Formosa, T. 
rot = Thalassiosira rotula, P. tric = Phaeodactylum tricornutum, H. triq = Heterocapsa trique- 
tra, P. min = Prorocentrum minimum, C. furc = Ceratium furcoides, G. rad = Golenkinia 
radiate, Chlor. sp. = Chlorella sp., S. cing = Staurastrum cingulum, A. ucr = Anabaena 
ucrainica, M. aer = Microcystis aeruginosa, A. flos = Aphanizomenon  flos-aquae, C. 
sperm = Cylindrospermopsis raciborskii 


growth rates at other temperatures, with the consequence that species replace each 
other (Fig. 2.4). 

This exponential temperature response of natural communities is usually 
expressed in terms of Qj, values or Activation energies E,, both rooted in chemical 


thermodynamics (van ‘t Hoff and Arrhenius equations). The temperature Qio is 
normally defined as 
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Fig. 2.4 Phytoplankton growth rate (d™') for a mixed community comprised of polar, temperate 
and tropical species. The mixed community response is based on Eppley (1972) 


Qio = (ee eta (2.5) 


HRef 


where ur and Href are the rate (e.g. growth) at temperature T and the reference 
temperature Tref (Celsius). Qio can be simplified to 


Qo = HRef +10 (2.6) 
HRef 


because it gives the rate increases for a 10 °C increase in T and is related to the 
parameter b of the exponential increase: Qjọ = e!%. Eppley’s curve thus 
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corresponds to a Qio of 1.88. Typical Qio values for biological processes are 
between 2 and 3. 
The Arrhenius equation is very similar and reads 


u= Aer? (2.7) 


where A is a pre-exponential factor (time), E, is the activation energy (J mol '), 
R is the universal gas constant (8.314 J mol ! K» and T is the absolute tem- 
perature (K). Sometimes the universal gas constant R is replaced by the Boltzman 
constant k (8.617 10° eV K`’) and then E, is expressed in eV (energy per mole- 
cule) rather than J mol '. For the temperature range of seawater, E, and Qj values 
are related via 


and Qig = el (™e)] (2.8) 


where T is again given in degrees Kelvin. Eppley’s Qio of 1.88 corresponds to 
activation energies of about 0.47 eV or 45 kJ mol! at 20 °C. One should realize 
that this is the optimal community temperature response, i.e. no other limiting 
factors. Apparent activation energies and Qjo values in the ocean are ~0.30 eV 
(29 kJ mol ') and ~ 1.5, respectively, close to that of Rubisco (Edwards et al. 
2016). 


2.2.3 Light 


Photosynthesis is a light dependent reaction, and light intensity has a major impact 
on growth rates. The relationship between photosynthesis and irradiance is nor- 
mally presented as a P versus E curve, where E refers to radiant energy (mol quanta 
m? s7’). Multiple equations have been presented to represent the photosynthesis to 
light relation, which differ in the number of parameters and whether or not they 
include the photo-inhibition effect at high light intensities or respiration of the 
autotroph. Photorespiration, the breakdown of photo-labile, intermediate carbon 
fixation products, is important in full-light exposed organisms, such as terrestrial 
plants, microphytobenthos and phytoplankton in the surface layer. 
Common simple limitation functions are the hyperbolic, Monod model: 


E 


fiim(E) = EFK) 


(2.9) 


where fim(E) is the light limitation function (0 < fiim(Œ)< 1) Kg is a 
light-saturation parameter (typically 50-150 pmol quanta m`? s™! for marine 
phytoplankton), and the Steele model (1962): 
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Siim(E) = E ¢(ha) (2.10) 


where Emax is typically 50-300 umol quanta m° s | for marine phytoplankton 
(Soetaert and Herman 2009). The Steele model represents both the initial increase 
and the subsequent decrease due to photo-inhibition with only one parameter 
(Fig. 2.5). 
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Fig. 2.5 Example of light inhibition functions in use. The Web and Monod models do not have 
light inhibition and only show different saturation behaviours, while the Steele and Platt models do 
incorporate the decrease in phytoplankton growth at high light levels, due to photo-inhibition 
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The Webb et al. (1974) model is based on an exponential: 
fim(E)=1-— E = eric) (2.11) 


where Pmax is the maximum rate at high light and o is the initial slope (increase in P 
with E at low light intensity). This equation ignores photo-inhibition. Alternatively, 
one can use the two-parameter Platt et al. (1980) equation: 


-BE 


fim(E) = 1 — [1 — lea) | [Fe] (2.12) 


where is p the intensity at the onset of photo-inhibition. Figure 2.5 illustrates the 
light limitation functions or PE curves presented above. 


2.2.4 Nutrient Limitation 


Growing phytoplankton needs a steady supply of resources to maintain growth. 
Nutrient uptake and growth kinetics are usually described using Monod or Droop 
kinetics. The former is the simpler model and normally used for steady-state 
conditions, while the Droop or internal quota model is preferred for transient 
conditions, e.g. in fluctuating environments. The equation for nutrient limitation 
following Monod kinetics is: 


S S 
H = Umax (S+K,) oF fiim(S) = (aa ) (2.13) 
where S is the substrate concentration of the medium water, fiim(S) is the nutrient 
limitation function, [max is the maximal growth rate, and K, is the half saturation 
constant for growth. 
The Droop equation expresses growth rate as a function of the cellular quota 
(Q) of the limiting nutrient (Droop 1970): 


4 Q 7 Omin 
H= hpa g (2.14) 
where Qmin is the minimum cellular quota for growth. Maximum growth rate on 
substrate (Hmax) and cellular quota (Hinar) are related Via ng, = La 


where Qinax is the maximum cellular quota if S increases. 
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2.3 From Theory and Axenic Mono-Cultures to Mixed 
Communities in the Field 


Progress in theory, creativity in experimental design, and dedicated hard laboratory 
work has generated process-based understanding of phytoplankton growth in the 
laboratory. This body of knowledge has deepened our understanding and guided 
our modelling efforts and field observation strategies, but we need to make many 
assumptions before we can apply this mechanistic approach to the field. 

Let us return to our master Eq. (2.1): P = u-B-fiim (resources, environmental 
conditions). Ignoring environmental conditions, such as temperature, and substi- 
tuting the simplest expressions introduced above we arrive at: 


E 
P= Hmax * B. j 2 (2.15) 
(E+Ke) (S+K,) 


This equation for primary production contains 6 terms that need to be quantified for 
the case of a single limiting nutrient and a single phytoplankton species. The light 
availability (E) and nutrient concentration (S) display spatial and temporal gradients 
in nature, and the maximum growth rate Umax and half-saturation dependences (Kg 
and K,,) require experimental or laboratory studies. 


2.3.1 Does Diversity Matter or Not? 


One of the most critical restrictions on the use of mechanistic complex models is 
related to phytoplankton diversity. Hutchinson (1961) identified the paradox that 
phytoplankton is highly diverse, despite the limited range of resources they com- 
pete for, in direct contrast to the competitive exclusion principle (Hardin 1960). 
Seawater typically contains tens of different species of primary producers, many for 
which there are no maximum growth data and known limitation functions. 
Accordingly, it is not feasible to simply apply Eq. 2.15 to individual species in the 
field and sum their contributions to obtain the primary production. Besides these 
theoretical arguments against the single species approach, there are also empirical 
reasons. Primary production and its dependence on environmental conditions 
(nutrients, temperature, light) are normally quantified at the community level in the 
absence of techniques to quantify species-specific primary production in natural 
waters. This discrepancy between, on the one hand, mechanistic, single-species 
approaches in the laboratory and, on the other hand, quantification of community 
responses and activities is somewhat unfortunate (Box 2.2). 
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2.3.2 Chi the Biomass Proxy 


The biomass of the primary producer (B) is the second term in our master equation 
and quantifying this term in natural systems is more difficult than one initially 
would anticipate. Particulate organic carbon (POC) concentrations (g C per unit 
volume) are a direct measure of phytoplankton biomass in laboratory settings with 
axenic cultures. However, in natural systems, the pool of particulate organic carbon 
comprises not only a mixture of phytoplankton species, each with its own maxi- 
mum growth rate, temperature, light and nutrient dependence, but also a variable 
and sometimes dominating contribution of detritus (dead organic matter), bacteria 
and other heterotrophic organisms. It is for this reason that chlorophyll concen- 
trations (Chl) are used as a proxy for living primary producer biomass. The 
rationale is that Chl is only produced by photosynthesizing organisms, degrades 
readily after death of the primary producers and can be measured relatively easily 
using a number of methods. Primary producer biomass (B) can then be calculated if 
one knows the C:Chi (or Chl:C) ratio of the phytoplankton. However, this ratio 
differs among species and depends on growth conditions, in particular light and 
nutrient availability (Cloern et al. 1995). Chl:C ratios vary from ~ 0.003 to ~ 0.055 
(gC gChlI'; Cloern et al. 1995), complicating going from phytoplankton growth to 
primary production. The very reason that Chl is such a good proxy for photosyn- 
thesizing organisms is also the reason why it is not well suited to the task of 
partitioning itself among different phytoplankton species: it is in all primary pro- 
ducers harvesting light energy. Accessory and minor pigments such as zeaxanthine 
and fucoxanthine, do, however, have some potential to resolve differences among 
phytoplankton groups, but not at the species level. 


2.3.3 Light Distribution 


The distribution and intensity of photosynthetically active radiation in seawater is 
governed by the intensity at the sea surface (Ep) and scattering and absorption of 
light, with the result that light attenuates with depth. The decline of light intensity E 
with water depth z can be described by a simple differential equation, expressing 
that a constant fraction of radiation is lost: 


dE 
— = —KparE (2.16) 
dz 


where the proportionally constant kpar is known as the extinction coefficient (m }). 
Solving this equation using the radiation at the seawater-air interface (Eg) yields the 
well-known Lambert—Beer equation: 


E = Ege” "rare (2.17) 
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The extinction coefficient kpar includes the absorption of radiation by water (kw), 
by the pigments from various primary producers (kcn), by coloured dissolved 
organic matter (kpoc), and by suspended particulate material (k,,,,). The light 
extinction coefficient of pure water (kw + 0.015-0.035 m~!) depends on the wave 
length of light, with longer wavelength (red) being adsorbed more strongly than 
shorter wavelengths (blue); this is the cause of the blue appearance of clear water. 
The other light extinction components have a different wavelength dependence: the 
attenuation coefficients of dissolved organic matter (kpoc; “gelbstoffe”’) and detritus 
(kspm) increase with shorter wave length, while that of phytoplankton (Kcy1) varies 
depending on the species, i.e. the pigment composition of the primary producers 
(Kirk 1992; Falkowski and Raven 1997). 

Oceanographers often divide ocean waters into two classes with respect to light 
absorption: case 1 waters in which phytoplankton (<0.2 mg Chl a m°) and its 
debris add only to kw, and case 2 waters which have high pigment concentration 
and light attenuation because of (terrestrially derived) dissolved organic carbon and 
suspended particulate waters. The overall light attenuation (kpar) in case 1 waters 
can be approximated by (Morel 1988): 


kpar = 0.121 x Chi*8 (2.18) 


where Chl is in mg Chl a m”. 


Other useful empirical relations link light attenuation (kpar) to the Secchi depth 
(Zsec, m), the depth at which a white disk disappears visually: 


Kpar = 3 ; (2.19) 
ZSec 


where q varies from 1.7 in case 1 waters to 1.4 in case 2 waters (Gattuso et al. 2006) 
and 


1.09 
ZSec 


kpar = 0.4 + (2.20) 


for turbid estuarine waters (Cole and Cloern 1987). 

Light attenuation coefficients vary from 0.02 m™' in oligotrophic waters, 
0.5 m™' in coastal waters, and to >2 m™' in turbid waters Light attenuation by 
water and phytoplankton dominate in the open ocean and on the shelf. In other 
coastal waters, including estuaries, phytoplankton and suspended particles domi- 
nate light attenuation, while light attenuation is primarily due to suspended particles 
in more turbid systems (Heip et al. 1995). 

The light attenuation governs the euphotic zone depth (Zeu, m), i.e., the depth 
where radiation is 1% of the incoming: 


4.6 
In 0.01 = —ŘPARZEU or ZEU = =~ (2.21) 
kpar 
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Fig. 2.6 Conceptual figure of euphotic zone depth (solid black line) as a function of suspended 
particulate matter (SPM, mg L~’) and phytoplankton concentrations (Chl, mg m °). The euphotic 
zone is more >200 m in the clearest ocean water with very low phytoplankton and light attenuation 
by water itself dominates. In most of the ocean, phytoplankton dominates light attenuation and 
euphotic zone depth scales with phytoplankton concentration. In estuaries and other turbid 
systems, dissolved organic matter and in particular suspended particles attenuate light and the 
euphotic zone narrows to less than one meter. Coastal systems and eutrophic parts of the ocean are 
in between. Case | and 2 oceanic waters are indicated. Light attenuation due to phytoplankton was 
modelled following Morel (1988; Eq. 2.18), while that due to suspended particles followed Cloern 
(1987) and euphotic depth was calculated as 4.6/kpar 


The euphotic zone is a key depth horizon in aquatic sciences because photosynthesis 
is largely limited to this zone. Moreover, the bottom of the euphotic zone is often used 
as reference for export of organic matter. Euphotic zone depths vary from about 
200 m in the oligotrophic ocean, to tens of meters in shelf systems, to meters in 
coastal waters and a few decimetres in turbid and/or eutrophic estuaries (Fig. 2.6). 


2.4 Factors Governing Primary Production 


Having presented the factors governing phytoplankton production in laboratory 
studies and the limitations in applying that knowledge to natural systems, we have 
all the ingredients to explore the factors governing the (depth) distribution and rate 
of primary production in natural ecosystems. 


2.4 Factors Governing Primary Production 23 
2.4.1 Depth Distribution of Primary Production 


Consider a system with a light profile following the Lambert—Beer equation (2.17) 
with E, = 10 mol md! and kpar = 0.1m! (corresponding to a euphotic zone 
of 46 m) and a nutrient pattern as shown in Fig. 2.7. Nutrients are low in the upper 
25 m (N = 0.1 pmol m °) and then exponentially increase with a depth coefficient 
0.1 m™' to a maximum of 10 umol m”. 

If we further assume (1) that physical mixing homogenizes phytoplankton 
biomass (B = constant), (2) that there is only one limiting nutrients (N), and (3) that 
light and nutrient limitations can be described by Monod relations with parameter 
Kg and Ky. This allows combining Umax and B into a depth independent maximal 
production P,,. The modelled P is then: 


E N 
P=!" EL Ke) NFR) or 


Taking Ky and Kg values of 1, i.e. 10% of E, and maximum N at depth, and 
combining Eq. 2.22 with the light and nutrient profiles, we can then calculate the 
primary production as a function of depth (Fig. 2.7, green curve). Although these 
light and nutrient profiles and the model parameters Ky and Kg numbers have been 
chosen arbitrarily, they are reasonable and generate a representative depth profile for 
primary production with a subsurface maximum, as observed as a deep chlorophyll 
maximum (Fig. 2.8). In the upper 25 m, primary production is rather low because of 
nutrient limitation and declines slightly with depth because of light attenuation 
(Fig. 2.7). Primary production is optimal at depths between 25 and 40 m, i.e. where 
the nutricline and the lower part of the euphotic zone overlap. Primary production 
below 25m is primarily light-limited, but accounts for about 75% of the 
depth-integrated primary production. Increasing surface-water nutrient concentra- 
tions or the phytoplankton affinity for nutrients (lowering Ky) would increase pri- 
mary production in the top 25 m, but not so much at depth (Fig. 2.9a). Increasing the 
photosynthetic performance at low light levels (lowering Kg) would increase pri- 
mary production at depth (Fig. 2.9b). Phytoplankton living in the surface ocean can 
thus optimize their performance by investing in nutrient acquisition, while those 
living in the subsurface would best optimize their light harvesting organs. This 
simple model explains why deep chlorophyll maxima occur in low-nutrient systems 
and why the depth distribution of primary production follows light in eutrophic 
systems (e.g. during early spring in Bermuda Atlantic station, Fig. 2.8). 


2.4.2 Depth-Integrated Production 


The overall control of light on depth-integrated production underlies 
satellite-derived algorithms for primary production and coastal predictive equations. 
For ecosystem and biogeochemical studies, the focus is on net primary production, 
i.e. carbon fixation minus phytoplankton respiration, expressed per m? and unit time 
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Fig. 2.7 Light and nutrient distribution versus depth and resulting primary production. The 
subsurface maximum often results in a deep chlorophyll maximum 


2.4 Factors Governing Primary Production 25 


Arbitrary Units 


0 0.2 0.4 0.6 0.8 

0 

50 

100 
G 
fa 
a 
[om 
1] 
(a 

150 

200 February 
August 
250 


9529 


Fig. 2.8 Fluorescence data at station BATS (Bermuda Atlantic Time Series) in February and 
August 1996 


(day, year). Behrenfeld and Falkowski (1997) showed that depth-integrated net 
primary production (P, g C m’ yr ') can be estimated as: 


P = Popi X Chl X Zeu X DL X fiim(E) (2.23) 
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Fig. 2.9 a The impact of nutrient availability on the vertical distribution of primary production. 
Low half-saturation constants Ky imply high availability of nutrient for phytoplankton. b The 
impact of light harvesting efficiency on primary production. A high affinity (low Kẹ) for light 
causes higher primary production at depth. The nominal run is presented in Fig. 2.7 


where Popt is the maximum daily photosynthesis rate (mg C (mg Chl) ' h7'), Zeu is 
the euphotic zone depth, DL is day length (h), and fiim(E) is a light limitation 
function. The similarity with our master Eq. (2.1) is evident, when nutrient limi- 
tation and environmental conditions are ignored. Integrating (Eq. 2.1) with depth to 
Zgu, and with time to sunset, we arrive at: 

P= SE Tonrise Hmax * B Sim(E) (22A) 
which is identical to (2.23), with Pop: = max; Chl =B, jj” = Zeu, and 
Isunrise = DL. 

Behrenfeld and Falkowski (1997) showed that 85% of the variance in global net 
primary production can be attributed to depth integrated biomass (Chl x Zeu) and the 
maximal photosynthesis parameter Pop, with other factors, such as differences in 
light limitation functions, depth distributions of phytoplankton biomass and day 
length (DL), being less important. Consequently, the most rudimentary model 
would be (Falkowski 1981): 


P=wx Chl x zgy x Eo (2.25) 


stating that net primary production (P) scales linearly with depth integrated biomass 
(Chl x Zgy), incoming radiation (E,) and an optimal photosynthetic parameter (W). 
Similar semi-empirical relations are often used in estuaries (Cole and Cloern 1987; 
Heip et al. 1995): 


P = a+b(Chl x zeu x Eo) (2.26) 
where a and b are regression coefficients that are system specific. 
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2.4.3 Critical Depths 


The overall governing role of light on primary production and phytoplankton 
dynamics also underlies the use of two critical depth horizons, often credited to 
Sverdrup (1953): the compensation depth (z,) and critical depth (Zer). These were 
introduced to understand and predict spring blooms in the ocean. At the compen- 
sation depth (Ze), phytoplankton photosynthesis is balanced by community respi- 
ration (Fig. 2.10), i.e. the depth of the radiation level at which photosynthesis (by 
phytoplankton) compensates their respiration (E,). This compensation depth should 
not be confused with the physics governed mixed-layer depth (Zma) and the critical 
depth (Zer), where primary production integrated through the water column and over 
the day will equal the daily water column integrated community losses of carbon 
(Sverdrup 1953; Fig. 2.10). These depths are pivotal to the formation of 
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Fig. 2.10 Production and respiration (left) and depth integrated production and respiration (right) 
as a function of water depth. The critical depth is where depth integrated production and 
respiration balance (200 m), the compensation depth is where production and respiration of 
phytoplankton balance (30 m). The euphotic zone is governed only by light attenuation and is 
46 m. Po is 30 mol m`? d !; Ro is 0.05*Po; Eg = 30 mol m7 d`!; kpag = 0.1 m'! 
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phytoplankton blooms in the oceans (Sverdrup 1953). If the mixed layer is deeper 
than the critical depth (z,,), then phytoplankton will spend relatively too much time 
at low irradiances and carbon losses are not compensated by sufficient growth. 
Conversely, if the mixed layer is shallower than Zer, phytoplankton communities 
can grow and blooms can develop. Assuming that carbon losses (R,) are constant 
with depth, there is no nutrient limitation, and gross primary production is linearly 
related to radiation, which in turn depends exponentially on depth (Eq. 2.18), 
primary production is described by: 

P = Poe *rk=, where Po is the surface productivity. One eventually arrives at 
following relations for Sverdrup’s critical depth, Zer: 


(1 = ekparzer ) Ec E Ro 


= = 2.27 
KPARZer Eo Po ( ) 


where FE. is the radiation level at the compensation depth and R, is the 
depth-independent community respiration rate (Sverdrup 1953; Siegel et al. 2002). 
Clearly, light attenuation is a major factor, not only governing Zeu, but also z, and 
Zer. The critical depth (Zer) is usually 4 to 7 times higher than the euphotic zone 
depth (Zeu). The compensation depth (z,) is typically 50-75% of the euphotic zone 
depth (Siegel et al. 2002; Sarmiento and Gruber 2006; Fig. 2.10). For simplicity, 
the compensation depth is often taken equal to the euphotic zone depth; this should 
be discouraged, because it implies that community respiration represents only 1% 
of maximal production. The depth of the euphotic zone (Zgy) is an optical depth 
governed by the light attenuation and thus only indirectly impacted by phyto- 
plankton via their effect on kpar, while the compensation depth depends on the 
community structure (algal physiology and heterotrophic community). The Sver- 
drup critical depth model is simple, instructive and predictive: it can explain bloom 
initiation when mixed layers shallow and link it to physical sensible and quantifi- 
able parameters. However, it is sometimes difficult to apply because of inconsis- 
tencies and uncertainties in the parameterisation (phytoplankton vs. community 
respiration and other phytoplankton losses) and the validity of the assumptions (no 
nutrient limitation, well-mixed layer). 

The critical depth horizon concept has been developed for deep waters, but a 
similar approach can be applied to shallow ecosystems. In shallow coastal systems, 
it is the relative importance of water depth and euphotic zone depth that governs 
(a) where production occurs and (b) whether phytoplankton biomass will increase 
or not. If water depth is less than the euphotic depth (Zguy) light reaches the seafloor 
and primary production by microbial photoautotrophs (microphytobenthos), as well 
as macroalgae and seagrasses, may occur. Gattuso et al. (2006) showed that this 
may happen over about 1/3 of the global coastal ocean. If water depth exceeds the 
euphotic zone by more than a factor 4-7 then phytoplankton losses in the dark 
cannot be compensated fully by photosynthesis and phytoplankton communities 
will lose biomass (Cloern 1987; Heip et al. 1995). Vice versa, if water depth <4—7 
times Zey phytoplankton growth is maintained. Consequently, shallowing of 
ecosystems (e.g. water flowing over a tidal flat or development of stratification) 
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Fig. 2.11 Conceptual 
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stimulates phytoplankton community growth, all other factors remaining equal, 
while deepening of water bodies will cause a decline. Moreover, in turbid systems 
where the light attenuation (kpar), and thus Zgy (Fig. 2.7), are governed by sus- 
pended particulate matter dynamics, phytoplankton communities may experience 
variable twilight conditions and have difficulty maintaining positive growth. Con- 
sequently, when turbid rivers and estuarine waters with high nutrients reach the sea, 
particles settle and light climate improves, phytoplankton blooms may develop and 
utilize the nutrients (Fig. 2.11). 

Sverdrup’s critical depth hypothesis is based on the assumption that phyto- 
plankton biomass and phytoplankton losses are homogenously distributed in the 
mixed layer. However, the mixed layer with uniform temperature as used in 
Sverdrup’s approach does not match with the layer of turbulent mixing in the ocean 
(Franks 2015). It is more realistic to represent phytoplankton biomass (B) as 
governed by the balance between production, respiration losses and transport by 
eddy diffusion and particle settling. Again, we assume gross primary production is 
linearly related to radiation (which declines exponentially); hence: P = Pye *™, 
Phytoplankton respiration loss is considered a first order process: Loss = rB with a 
first-order rate constant (r). Under the assumption of steady-state we then arrive at 
(see Box 1.1): 


@B dB eae 
K Z = Wiz = rB = Poe PARS: (2.28) 
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where K, is the vertical eddy diffusion coefficient (m’s '), w is the settling velocity 
(m s}; positive downwards), the other terms have been defined before. Consid- 
ering a semi-infinite domain, i.e. B = 0 at large depth, and phytoplankton biomass 


Bo at the water-air interface, we obtain the following solution: 


Po az Po —kparz 
B = Bo 7 E + 7 e PAR (2.29) 
K-kpar + wkpar =F K-kpar + wkpar =i 


. w—4/w? +4rK; 


The second exponential term accounts for light-dependent production, while the 
first exponential comprises water-column mixing, phytoplankton settling, and 
phytoplankton losses. To simplify matters, we assume that phytoplankton biomass 
is zero at the air-water interface. The first and second term then balance if 


w—/ w? +4rK; 


= —kpar. After re-arrangement to isolate the eddy diffusion coefficient, 


2K, 
we obtain 
—k 
K = (2.30) 
kar 


In other words, the vertical eddy diffusion coefficient K, should be less than “pan 
PAR 


for positive values of phytoplankton biomass (B). 

Huisman et al. (1999) presented a more elaborate model on phytoplankton 
growth in a turbulent environment, including a feedback between phytoplankton 
biomass and kpar. Through scaling and numerical analysis of a model without 
phytoplankton sinking (w = 0), they derived a relationship between the maximum 


turbulent mixing coefficient K, and Kpar: K; = oat If we also ignore phyto- 
PAR 
plankton advection (w = 0 in Eq. 2.30), K,< g—, fully consistent with Huisman 
PAR 


et al. (1999). The critical turbulence level for phytoplankton growth is thus 
inversely related to the square of the attenuation of light. Moreover, the phyto- 
plankton loss is the scaling factor. For turbid systems such as estuaries and other 
coastal systems with high light attenuation (kpar), turbulent mixing should be 
minimal to allow net growth, consistent with observations by Cloern (1991) that 
phytoplankton blooms develop during neap tide when turbulent mixing intensity is 
lowest. Conversely, in clear, oligotrophic waters, light attenuation is limited and 
phytoplankton blooms can occur at relatively high mixing rates. Sinking phyto- 
plankton (w > 0) will lower the numerator of Eq. 2.30 and thus lower the critical 
turbulence levels, while buoyant phytoplankton (w < 0) will increase the maximal 
allowable turbulence, and thus the scope for phytoplankton growth. 
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Fig. 2.12 Production as a function of individual biomass for a wide variety of organisms, 
including animals, plants and prokaryotes (based on data compiled by Hatton et al. (2015)) 


Box 2.1: Phytoplankton size based traits 

The intrinsic maximum growth of phytoplankton varies with size (Fig. 2.2). 
Metabolic activity of organisms usually scales with size and when expressed 
in terms of mass or volume (V) follows a simple power law 


u= aV”, (2.31) 


where b = —0.25 according to the metabolic theory of ecology (Brown et al. 
2004). Accordingly, the smaller the organism, the higher the intrinsic max- 
imum growth rate. This power law relationship holds over orders of mag- 
nitude and across a wide range or organisms (autotroph and heterotroph, 
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eukaryotes and prokaryotes; e.g., Fenchel 1973) and implies that smaller 
organisms have the highest intrinsic growth (Fig. 2.12). 

However, some cell components are non-scaleable, such as the genome 
and membrane, and consequently this power-law appears to break down in 
the range of nanoplankton (2-20 um). There is a trade-off between the size 
dependence of physiological traits (Ward et al. 2017). Burmaster’s (1979) 
equation can be used to illustrate this: 


Hmax * Dsize 
ye ee 2.32 
ene Umax i Onn F Ore ( ) 


where the maximum growth for a certain size (Msize) depends on maximum 
nutrient uptake (Osize), minimum cell quota (Qmin) and theoretical maximum 
growth rate (Umax). Maximum nutrient uptake and requirement per cell scale 
positively with cell size (Fig. 2.2b, dashed blue line), while theoretical 
maximum growth rates scale negatively (Fig. 2.2b solid blue line). The result 
is an optimum in growth rate for phytoplankton in the nanoplankton range 
(Fig. 2.2b, black line). Very small picoplankton cells have a low intrinsic 
growth rate that will increase with size because more volume is then available 
for catalysing and synthesizing. The intrinsic growth rate of microplankton 
cells will decrease with increasing size, as with most organisms, for multiple 
reasons, including the increase in intracellular transport distances between 
cellular machineries (Marañón et al. 2013). 


Box 2.2: Phytoplankton diversity, rate measurements and biogeo- 
chemical models 

The high number of different species in each water sample poses a challenge to 
link the species-specific growth parameters obtained in the laboratory with 
measurements of phytoplankton growth in the field and modelling of phyto- 
plankton primary production for natural, mixed communities. Gross primary 
production is normally quantified by the production of oxygen, using either 
'80-labelling or the differential evolution of oxygen in light and dark. The most 
common technique for measuring primary production is the BG labelling tech- 
nique, but this method provides a result in between gross and net photosynthesis, 
depending on the duration of the incubation. Both approaches quantify primary 
production for the total community, rather than for specific species. Biological 
oceanographers have developed methods to quantify group-specific primary 
production, based on dilution approaches or the incorporation of isotopically 
labelled bicarbonate into biomarker or flow-cytometry separated groups of 
organisms (Laws 2013). These group-specific primary production measurements 
can be compared more directly to laboratory data. 


2.4 Factors Governing Primary Production 33 


Biogeochemical modellers have explored a number of strategies to 
incorporate differences among phytoplankton species into their ecosystem 
models; e.g. the plankton functional group approach and phytoplankton size 
or trait based approaches. The former approach is limited to a few plankton 
groups that are representative for certain biogeochemical fluxes (e.g. N2- 
fixers, diatoms, small and large phytoplankton, coccoliths; Sarmiento and 
Gruber 2006). The size-based approach makes use of the systematic rela- 
tionships between phytoplankton size and activity (e.g. Fig. 2.2), but some 
processes do not scale in a simple way with size. Trait- and genome-based 
approaches are the most recent, and they consider emergent phenomena 
(Follows et al. 2007). These approaches are instructive and needed to further 
our understanding and predictive capabilities in times of global change, but 
they are so far difficult to link with observations in the field. 
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The Return from Organic to Inorganic 3 
Carbon 


Almost all the organic matter produced is eventually consumed and respired to 
inorganic carbon because organic matter preservation via burial in accumulating 
sediments (~ 0.2-0.4 Pg y_') represents only a very small fraction of that produced. 
Global phytoplankton production is about ~50 Pg C y ', while phytoplankton 
biomass is ~1 Pg, implying a turnover of one week (0.02 y_'). Marine macro- 
phytes have a similar global biomass, but a production of only 1 Pg C y |; the 
turnover is thus ~1 y ' (Smith 1981). These high turnover rates (compared with 
global terrestrial vegetation turnover of about one to two decades, Field et al. 1998) 
imply not only steady production, but also efficient consumption of organic matter. 
There are multiple organic matter loss pathways (respiration by autotrophs and 
heterotrophs, grazing, viral lysis, detrital route), but all eventually result in respi- 
ration and release of inorganic carbon. 

Biogeochemists constructing carbon budgets normally lump together the various 
organic matter loss pathways and focus instead on the quantification of the organic 
carbon to inorganic carbon transformation. However, for a more detailed under- 
standing, the elucidation of the link with other biogeochemical cycles and the 
identity of organisms involved, the various pathways have to be resolved. Before 
discussing the processes and mechanisms involved in these routes, it is instructive 
to distinguish between living organic matter that has the capability to reproduce 
(primary and secondary producers) and dead organic matter (i.e. detritus). While 
living organic matter is, by definition, fresh and thus labile, detrital organic matter 
pools represent a heterogeneous mixture of compounds from various sources which 
have been mixed together and which may have distinct compositions, degradation 
histories and reactivities (Box 3.1). Another distinction is between organic matter in 
the dissolved and particulate phases. Concentrations of DOC (order 10- 
100 mmol m™*) are usually one order of magnitude higher than those of POC 
(order 1-10 mmol m °) in the euphotic zone of the ocean, while sediment POC 
(order 10° mmol m™°) concentrations are normally three orders of magnitude higher 
than those of DOC in pore water (100-1000 mmol m™°). Moreover, the DOC pool 
is detrital (with viruses), while the particulate organic carbon pool usually 
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represents a mixture of living organisms, their remains, and other detrital organic 
carbon inputs. As we discussed in Chap. 2, the POC concentrations in the water 
column can be a poor proxy for phytoplankton carbon. This distinction between 
dissolved and particulate pools is pivotal for understanding the fate of organic 
matter. Dissolved organic matter is transported as a solute with the water, while 
particulate organic matter is subject to gravity, which results in settling of organic 
particles. Consequently, the residence time of DOC in an aquatic system is much 
longer than that of POC. Consumption of dissolved organic matter primarily 
involves microbes and other small organisms, which use it for energy and nutrient 
acquisition, while larger organisms generally prefer particulate organic matter. 
Microbes use extracellular enzymes to solubilize particulate organic matter before 
they can utilize it. 

In this chapter, we discuss carbon consumption in the context of the “biological 
pump”. The latter depends on the fraction of primary produced organic matter that 
survives degradation in the euphotic zone and that is exported from surface water to 
the ocean interior, where it mineralized to inorganic carbon, with the result that 
carbon is transported against the gradient of dissolved inorganic carbon from the 
surface to the deep ocean. This transfer occurs through physical mixing and 
transport of dissolved and particulate organic carbon, vertical migrations of 
organisms (zooplankton, fish) and through gravitational settling of particulate 
organic carbon (Volk and Hoffert 1985; Sarmiento and Gruber 2006). We first 
discuss carbon consumption pathways in the euphotic zone, then the factors gov- 
erning export of organic carbon from the euphotic zone and particle degradation in 
the ocean interior. Processes specifically related to organic carbon consumption in 
the coastal zone are presented in Box 3.2. 


3.1 Carbon Consumption Pathway in the Euphotic Zone 


For pelagic ecosystems, Legendre and Rassoulzadagan (1995) proposed a contin- 
uum of trophic pathways with the herbivorous food-chain and microbial loop as 
food-web end members (Fig. 3.1). The classical linear food-chain end-member 
involves grazing by zooplankton on larger phytoplankton and subsequent predation 
on zooplankton by either larger zooplankton or another predator. In such a linear 
food-chain a predator can either lead to high phytoplankton biomass (in a system 
with phytoplankton, herbivore and a predator) or reduced phytoplankton biomass 
(in a system with four levels). Changes in predator abundance can, thus, lead to 
trophic cascades (Pace et al. 1999). The microbial loop end-member involves not 
only phytoplankton, as basal resource, but also dissolved organic carbon (Azam 
et al. 1983). Dissolved organic carbon is used by heterotrophic bacteria for growth 
and respiration, and these heterotrophic bacteria are, in turn, consumed by micro- 
zooplankton (20-200 um; ciliates, radiolarian, foraminifera) that are predated upon 
by larger zooplankton. Consequently, dissolved organic carbon is transformed, via 
a bacterial-microzooplankton loop, to zooplankton. These two end-member carbon 
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processing pathways are connected at multiple levels. Small phytoplankton can be 
consumed directly by microzooplankton. Dissolved organic carbon is produced in 
multiple ways and by various organisms, both by primary producers and consumers 
of organic carbon (Fig. 3.1). DOC release by primary producers occurs passively 
by leakage and actively during unbalanced growth during nutrient limitation 
(Anderson and LeB Williams 1998; Van den Meersche et al. 2004). Another direct 
pathway from phytoplankton to dissolved organic pool involves viral lysis (Suttle 
2005). Viruses are a major cause of phytoplankton mortality in the ocean, partic- 
ularly in warmer, low-latitude waters. Sloppy feeding by herbivores and incomplete 
digestion of prey by consumers are other sources of dissolved organic carbon. 
Heterotrophic microbes use extracellular enzymes to solubilize particulate organic 
carbon and use this and other dissolved organic carbon resources for growth and 
maintenance. Part of the microbial heterotrophic production is used by microzoo- 
plankton; another part of the heterotrophic community is subject to intense viral 
lysis and this causes release of dissolved organic carbon again. The efficiency of the 
microbial loop depends on multiple factors but in particular on the relative 
importance of predation and viral lysis to the mortality of heterotrophic microbes. 
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Fig. 3.1 Conceptual diagram of food web structure in euphotic zone. The linear food chain large 
phytoplankton-herbivore-predator (on the left with red arrow connections) has fewer levels than 
one with small phytoplankton at the base. The microbial loop refers to the flow from the dissolved 
organic carbon (DOC) via heterotrophic bacteria (Het. Bac.) and microzooplankton to predatory 
zooplankton (on the right with black solid arrows). Viruses play a major role in the mortality of 
phytoplankton and heterotrophic bacteria, and recycle organic carbon back to the DOC pool. Other 
sources of dissolved organic carbon (also dashed black arrows) includes exudation, sloppy 
feeding, etc. Particulate detritus pools and fluxes are not shown for simplicity 
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3.2 Factors Governing Export of Organic Matter 


While net primary production (gross production minus respiration by the auto- 
trophs) as presented in Chap. 2 is the property of interest for structure and func- 
tioning of food webs in the surface ocean, new and export production are most 
relevant for the functioning of deep-sea and sedimentary ecosystems and the role of 
the ocean in the global carbon cycle (Fig. 3.2). New production, also known as net 
community or ecosystem production, refers to net primary production minus the 
consumption of organic carbon by heterotrophs in the euphotic zone. At steady 
state, this new production should be the same as what is exported, i.e. export 
production. During periods of phytoplankton blooms new production may tem- 
porarily exceed export out of the euphotic zone, while rates of export are tem- 
porarily higher than new production during the senescence of phytoplankton 
blooms. 

The relations between export, new and net primary production are often 
expressed in ratios. The e-ratio is defined as: 


. export production 
e—ratio = 


(3.1) 


net primary production 


Export P 


Fig. 3.2 The various types of productivity in the ocean. GPP (gross primary production) is the 
total rate of carbon fixation/oxygen release by phytoplankton. Part of the GPP is used for 
autotrophic respiration (Ra) by the phytoplankton community, the other part represents net 
primary production (NPP). The organic matter produced is consumed by heterotrophs for growth 
and respiration (Ry) and the remaining organic matter is new and available for export (NEP, new 
production, export production). The net primary production (NPP) is most relevant for euphotic 
zone food web functioning, while the export production is relevant for the biological pump and 
fuelling deep-sea and benthic food webs. The export or e-ratio divides export/new production by 
net primary production (NEP/NPP) 
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where export production is often quantified via the use of sediment traps (Suess 
1980). Note that e-ratios based on sediment trap might be biased if dissolved 
organic matter is transported downwards by advection and/or eddy-diffusion and if 
zooplankton and other mobile vertical migrating organism contribute significantly 
to the total downward flux of organic carbon. An alternative method to quantify the 
relative importance of new to total production is the f-ratio, based on nitrogen use 
and recycling: 


new production 


(3.2) 


ratio = 
f new production + recycled production 


where new production is assumed to equal nitrate uptake and recycled production is 
taken equal to ammonium uptake as determined from '°N tracer experiments 
(Dugdale and Goering 1967). The rationale is that nitrogen availability limits pri- 
mary production, that nitrate is supplied to the euphotic zone from deeper waters 
(thus new to euphotic zone) and that recycling of organic nitrogen within the 
euphotic zone generates ammonium. Export of organic nitrogen should eventually, 
on the global scale, be equivalent to nitrate supply from the deep to surface ocean. 
However, nitrification (ammonium oxidation to nitrite and nitrate) in the euphotic 
zone, atmospheric deposition, nitrogen fixation and other nitrogen cycling pro- 
cesses (bacterial uptake of dissolved inorganic nitrogen, dissolved organic nitrogen 
generation and use) complicate this simple depiction. At steady state e and f-ratios 
should be the same if the underlying assumptions are correct. 

These export ratios vary from ~0.04 to ~0.72. Although multiple factors 
contribute to this range in export ratios (phytoplankton size, community structure, 
e.g. see Fig. 3.1), temperature and total primary production together account for 
87% of the variance in export production and ratios (Laws et al. 2000, 2011; Dunne 
et al. 2005). The empirical model of Laws et al. (2011): 


0.43T 
e—ratio = 0.04756 (078 = n pet (3.3) 


shows that e-ratios depend negatively on temperature (T, in Celsius) and scale with 
the power ~0.3 to net primary production (P). Export production consequently 
scales to net primary production with a power of ~ 1.3. Temperature is the single 
most important factor (Fig. 3.3), because organic carbon consumption by hetero- 
trophs is more temperature sensitive than light and nutrient-limited primary pro- 
duction by phytoplankton (Laws et al. 2000). The impact of primary production on 
the export ratio can be explained by the increase in the size (and thus settling rate) 
of phytoplankton with increasing primary production (Dunne et al. 2005). The 
global average export ratio is ~0.2 with eutrophic, high-latitude systems having 
high e-ratios and oligotrophic, low-latitude systems having the lowest ratios. 
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3.3 Particulate Organic Carbon Fluxes in Ocean Interior 


The efficiency of the biological pump depends not only on the rate of primary 
production (Chap. 2), the efficiency of carbon consumption within and export out 
of the photic zone (Fig. 3.3), but also on the depth at which the organic carbon is 
respired, because this determines the period during which carbon will be removed 
from the atmosphere (Yamanaka and Tajika 1996). 

This pivotal role of carbon transfer to and within the ocean interior has stimu- 
lated observational programs to measure particle fluxes using sediment traps. Suess 
(1980), Pace et al. (1987) and Martin et al. (1987) were among the first to explore 
such sediment trap data for the global ocean, and they observed a steady decline of 
organic carbon fluxes with depth because of degradation during settling (Fig. 3.4a). 
The empirical open ocean composite from Martin et al. (1987) is most often used: 


F, = F, (=) ý (3.4) 


Zo 


where z is water depth (m), Zo is the reference level for export (e.g., bottom of the 
euphotic zone), F,/,, are the fluxes of organic carbon (mol m`? d7’) at depths z and 
Zo, and b is a fitted coefficient with a value of 0.858. Although we have made much 
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Fig. 3.3 Relation between e-ratio (export production over net primary production) as a function 
of water temperature (Data compiled by Laws et al. 2000) 
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Fig. 3.4 a Fluxes of particulate organic carbon in the ocean interior based on the Martin et al. 
(1987), Suess (1980) and Pace et al. (1987) relations and data compiled by Lutz et al. (2002). Note 
that F,, has different meanings for the original equations but has been normalized to 1 for this plot. 
b Fluxes of particulate organic carbon in the ocean interior based on the Martin et al. (1987), Suess 
(1980) and Pace et al. (1987) relations and exponential curves with attenuation coefficient of 0.01, 
0.001 and 0.0001 (m™~). Martin et al.: F, = F», (2) CEE oK Fi 
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progress in our understanding of the governing processes, the basic early findings 
have not been falsified. The open ocean composite from Martin et al. (1987) is used 
in many earth system and ocean biogeochemical models. The parameter b typically 
ranges from 0.6 to 1.3, with lower values in low oxygen settings and higher values 
in productive high-latitude systems (Berelson 2001). This variance in b values has 
been attributed to community structure, oxygen, temperature and mineral ballasting. 
Alternative more complex formulations have been derived. 

To understand the factors governing particulate organic matter fluxes in the 
ocean interior, we will introduce a simple model for POC fluxes. At any water 
depth (z, directed downwards) below the export reference level (zo), the particulate 
organic carbon concentration is governed by the balance of particulate organic 
matter supply by settling particles and degradation: 


dPOC 
z 


w = —kPOC (3.5) 
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where w is the settling velocity (m d-'), POC is the concentration of particulate 
organic carbon (mol m °) and k is a first-order rate constant (d '). Assuming a 
constant first-order rate constant k (i.e. one-G model; Box 3.1), a constant settling 
rate (w), and a known export flux of organic carbon at the reference depth (F,, ) the 
solution is an exponential: 


Fy, sits 
POC, = Sew ea) (3.6) 


Moreover, F, = POC,w because w is constant, and thus, the flux of particulate 
organic carbon is given by: 


F, = Fe 0) (3.7) 


Accordingly, vertical organic carbon fluxes in the ocean are expected to decline 
exponentially with water depth and the attenuation is governed by the ratio between 
the organic matter degradation rate constant k and particle settling velocity w. 

Figure 3.4b presents the Suess, Pace and Martin curves again but now together 
with three curves based on the exponential mechanistic model with k/w values 
varying from high attenuation (0.01) to very low attenuation (0.0001). While none 
of these can describe the observations, the former agrees with the upper part, while 
the latter agrees with the lower part. This was already noted by Martin et al. (1987) 
and implies that the assumption of constant k and/or w is likely not correct. 

If our assumption of constant k is the main and only reason, then we would 
expect much improvement in model prediction with a multi-G or reactive contin- 
uum model (see Box 3.1). Figure 3.5 presents model predictions for both a multi-G 
and reactive continuum model based on model parameters for the phytoplankton 
decay experiments of Westrich and Berner (1984), presented by Boudreau and 
Ruddick (1991). While accounting for a decreasing reactivity of organic matter with 
depth improves the model performance in the upper part, there is still a large 
discrepancy at depth. However, better agreement can be obtained by fitting rather 
than imposing the organic matter degradation rate parameters based on literature 
values. The resulting gamma shape parameter (v > 0.6, depending on sinking 
velocity w chosen) is within the range of reported values (Arndt et al. 2013). Such a 
high v implies that settling organic matter is dominated by labile fractions, which 
seems reasonable. 

Alternatively, the particle sinking velocity may be a function of water depth. 
Particle sinking velocity (w) is governed by Stokes’ law: 


B 2gr°Ap 
w= On (3.8) 


where g is the gravitational acceleration (9.81 m? s7!), r is the particle radius (m), 
Ap is the density difference between seawater and the particle (kg m™°) and y is the 
dynamic viscosity of water (N s m ”). In other words, particle settling is governed 
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Fig. 3.5 Fluxes of particulate organic carbon in the ocean interior based on the Martin curve and 
model predictions using different parameterisation for organic matter degradation. The one-G has a 
fixed k of 0.03 (q™’», the three G-model and reactive continuum (RC) parameters were taken from 
Boudreau and Ruddick (1991), the two G-model has similar k parameterisation as the three-G 
model but without a refractory fraction. The RC fit curve was based on fitting parameter v to the 
Martin curve (v = 0.59). Such a high v implies a dominance of fresh material 


by the balance between gravitational acceleration of the particle (g x Ap) and the 
drag acting upon it via surface area (1°) and the friction of the fluid (u). Particle size 
and density are the crucial factors (Sarmiento and Gruber 2006). The density dif- 
ference Ap is much larger for calcite and clay minerals with densities of ~ 2700 
kg m ° than for organic matter with a density of ~ 1060 kg m °, a little higher 
than that of water ~ 1027 kg m°. An organic particle with a diameter of 100 um 
will settle ~12m d', while a pure calcite or clay mineral particles will 


sink ~600 md’. Decreasing the size of the particle by a factor 10 lowers sinking 
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velocities by a factor 100, because of the ir dependence, indicating that small 
reactive organic particles will never reach the seafloor. 

Although Stokes’ law is based on physical principles and generally valid, 
application to particles in the ocean is not straightforward. Particles are diverse and 
made up of aggregates of various plankton types, their detritus, faecal pellets and 
different mineral phases with distinct densities. Stokes’ law as formulated above is 
for spherical particles and most marine particles, organism and their remains may 
deviate strongly from this shape, with the consequence that they are somewhat less 
dependent on size. Moreover, particles sizes in the ocean vary over orders of 
magnitude. Size spectra analysis show that most mass occurs in the smaller size 
classes and cannot directly be linked to that of primary producers size spectra. 
During settling, particles interact, disaggregate and coagulate, with the consequence 
that the number and size distribution of particles changes. Natural organic matter 
aggregates have high porosity, which decreases the density difference between 
water and particles and thus gravitational acceleration. Nevertheless, settling 
velocities of natural particles have shown to follow a power law with particle size 
(Clegg and Whitfield 1990; Sarmiento and Gruber 2006). In situ settling velocities 
have been quantified by high-resolution underwater video systems and inferred 
from the time-lag of pulse arrival between sediment traps at different water depth in 
the same area. The latter approach showed that particle settling velocities increased 
with depth (Berelson 2002). 

To explore the impact of increasing settling velocities a linear and exponential 
increase of particle sinking rates with water depth was assumed (Fig. 3.6b) and 
Eq. 3.5 solved for these situations (Fig. 3.6a). A linear increasing rate of particle 
settling from 10 m d™' at the bottom of the euphotic zone to about 185 m, com- 
bined with a first-order degradation rate constant (k = 0.03 d; ~11 y~’) agrees 
with the parameterisation of Suess (1980) and Martin et al. (1987). Evidently the 
initial assumption of constant first-order rate constants and sinking velocities might 
be questionable. 

These agreements between the Martin’s ocean composite model and the reactive 
continuum (with high v) and linear increasing settling rates are not unexpected: it is 
in the equations. The flux equation for the reactive continuum model with uniform 


w: 
aw \* 

F, =F 3.9 

y (=) ( ) 


approaches the Martin et al. curve if aw << z. Similarly, the flux equation for 
constant reactivity but with linear increasing velocity according to wz = c z (c in 
d7’) is: 


F, = Fo E (3.10) 
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Fig. 3.6 a Fluxes of particulate organic carbon in the ocean interior based on the Martin curve 
and model predictions using different parameterisation for particle sinking rates. b Parameterisation 


for particle sinking rates in the ocean used to generate panel a. The solution to Eq. 3.5 with an 
-k 


å : é š re ‘ xk z +, @-b(z-z0) \ vob i 
exponentially increasing sink rate w = vo + vie ™ is: F, = F,,e% (e=z0) (ee) and with 


=k 


a linear increasing sink rate w = c z, it is: F, = F,, (4) j 


which returns Martin’s power law if = is —0.858. With a k value of 0.03 d! (11 
y`’), settling velocities should then increase with a c value of 0.035 from 10 to 
about 185 m d™' in the deep ocean, within the range of observed values (Berelson 
2001). 

Most of our knowledge of the nature, dynamics and magnitude of vertical 
organic particle fluxes in the ocean interior is based on sediment traps, which are 
one-side open vertical cylinders or conical containers that are fixed to a mooring or 
floating. These data may be biased because of hydrodynamics, degradation of 
material in the trap before recovery, and collection of swimming, foraging animals 
that are killed while feeding on the poisoned traps (Sarmiento and Gruber 2006). 

An alternative approach to quantify the depth attenuation of particulate organic 
carbon fluxes in the ocean involves the use of sediments as the ultimate sediment 
trap. The idea is that eventually all organic material arriving at the seafloor is 
respired and can be measured by sediment oxygen uptake (the rationale for this 
approach will be explained in Chap. 4). Our mass-balance approach is shown in 
Fig. 3.7. Assuming that sediment respiration at a particular depth reflects the par- 
ticulate organic carbon flux at that depth at any place in the ocean, we can calculate 
the total organic carbon flux as a function of depth by using the total area of the 
ocean at that depth. In other words, we assume a lateral homogenous ocean with no 
gradients in export between coastal and open ocean systems and between 
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Fig. 3.7 Schematic picture for the carbon balance for a depth layer in the interior ocean. The 
incoming flux of settling carbon (Fin) over the total surface area (Aot) of the ocean is balanced by 
the flux of carbon passing by (Fout) for the total surface area minus that of sediments (Asea) and the 
carbon that is respired in the water column (R,,) and in the sediment (R,). Moreover, sediment 
carbon respiration equals incoming flux (Fin) times global sediment area (Asea) at that depth 


high-latitude and low-latitude systems. Sediment oxygen consumption data of 
Andersson et al. (2004) were fitted with a power law and combined with ocean 
hypsometry (ETOP1) data. Moreover, by differencing carbon fluxes between two 
depth layers, we obtain total respiration in that depth zone. By substracting sedi- 
ment respiration from total respiration we can derive the water-column respiration. 
A similar approach can be applied to the Suess and Martin relationships, as already 
shown by these authors. Water-column respiration rates estimated from sediment 
oxygen consumption are consistent with respiration rates calculated from the Martin 
et al. (1987) and Suess (1980) relations (Fig. 3.8), both in terms of profile shape, as 
well as magnitude. These are completely independent estimates: the sediment 
approach is based on sediment oxygen consumption, while the Martin and Suess 
approaches are based on primary production and sediment trap data. Depth inte- 
gration of these global respiration versus depth profiles allows calculation of 
water-column respiration in the ocean interior (Table 3.1). Ocean interior 
water-column respiration estimates inferred from sediment oxygen consumption are 
intermediate between those derived from the Martin and Suess relationships (the 
latter is based on a global primary production of 50 Pg C y-'), and consistent with 
literature estimates (Dunne et al. 2007). However, part of the settling carbon flux 
arrives at the ocean floor and this is not included in the Martin and Suess relation- 
ship. Total respiration in the ocean interior is, therefore, substantially higher 
(Table 3.1). While water-column respiration declines systematically with increasing 
water depth, global sediment respiration rates increase again below 1—2 km because 
of ocean hypsometry: large parts of the ocean have water depths between 3 and 
6 km. As a consequence, at water depths more than 3 km sediments dominate 
organic carbon degradation, and thus oxygen consumption and carbon dioxide 
production (Fig. 3.8) This trend is robust because other sediment oxygen 
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Fig. 3.8 a Water-column respiration in the ocean interior calculated from the Suess and Martin 
relationships and from sediment oxygen consumption data and the mass balance model of Fig. 3.7. 
10 b Depth distribution of total respiration, water-column and sediment respiration based on a 
power fit to the Andersson et al. (2004) data. Sediment respiration dominates at water depths more 


than 3000 m because of hypsometry of the ocean. The Martin relations for respiration: Ry = 


13.1). The Suess flux relation F, = oR CO was combined with a global net 


primary production of 50 Pg C y '. Sediment data are from Andersson et al. (2004) and were fit 
with a power law: F, = 84.88(z)~°** (mmol C m? q» 


Table 3.1 Respiration in the ocean interior (Pg C y_') 


Water Water Water Sediment Total Total 
column column column Rs° respiration? Glud 
Martin* Suess? Rw| dataf 
100m 5.85 17.40 8.79 1.77 10.56 14.07 
200 m 3.05 8.64 5.53 1.59 WA 8.26 
1000 m 0.55 1.35 1.48 1.38 2.86 2.40 
2000 m 0.19 0.45 0.60 1.26 1.86 1.36 


“Water column respiration based on Martin et al. (1987) 

Water column respiration based on Suess (1980) and global primary production of 50 Pg C y~! 
“Total, sediment and water column respiration based on Andersson et al. (2004) data and approach 
outlined in Fig. 3.8 


“Total respiration based on Glud (2008) sediment data and same as approach as for Andersson 
et al. (2004) data 


consumption versus water depth curves generate similarly shaped curves and benthic 
contribution to deep-sea respiration (Table 3.1). Clearly, a more detailed look at 
sediment carbon processing is warranted, and this is provided in the Chap. 4. 
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Box 3.1: Organic matter reactivity 

Detrital organic matter represents a very heterogeneous pool of thousands of 
organic compounds. Individual compounds may originate from various 
organisms, come from different environments, may have been produced 
recently or thousands of years ago, and have experienced a different history 
before eventually making up the compound pool of dissolved or particulate 
organic matter. Moreover, for particulate organic matter these compounds 
may be associated with different mineral phases or incorporated in various 
ways into minerals (Arndt et al. 2013). All these factors and intrinsic dif- 
ferences in reactivity among organic compounds result in large differences in 
organic matter reactivity (Middelburg 1989). Organic matter reactivity is 
usually expressed in terms of a first-order rate constant (k), which originates 
from the one-G model of marine organic matter degradation (Berner 1964): 


dG», 
dt 


= —kGm, (3.11) 


where Gm is the concentration of metabolizable organic matter, t is time and k 
is a first-order rate constant (time ') assumed to be constant in time. This 
equation implies that the concentration (G,,) and rate Iu decrease expo- 
nentially with time. This model has been used successfully in various envi- 
ronments, from soils to sediments, from algal decomposition experiments to 
sewage degradation. However, there are two issues with its use (Middelburg 
1989). One, it assumes that we know a priori the partitioning between 
refractory organic matter (G,) and degradable organic matter (Gm) that add up 
to the total concentration (G). Two, organisms have been shown to prefer- 
entially utilize organic substrates, the more reactive being consumed first. 
This has a few consequences: the reactivity of the remaining organic matter 
decreases with reaction progress, i.e. time, and there are systematic changes 
in the composition of organic matter (Dauwe et al. 1999; Chap. 6). 

To account for the observed decline in reactivity of organic matter with 
time, two alternative approaches have been presented: the multi-G and 
reactive continuum models. The multi-G model (Jørgensen 1978; Westrich 
and Berner 1984) divides the organic matter pool into a discrete number of 
fractions (usually two reactive and one refractory) with a different reactivity, 
each of which undergoes first-order decay. The relevant equations are: 


G=) Gi (3.12) 


dG; 
dt 


= kG; (3.13) 
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dG 
—=-Y) kG; 14 
7 > G (3.14) 


where G; is the concentration of organic carbon in each group i, k; is the 
first-order reactivity of each group and i > 1. The one-G model (Eq. 3.11) is 
an example where i = 2, one labile class and one refractory group (with a k; 
of zero). Selective removal according to the reactivity of each group accounts 
for the decrease in reactivity, amount and rate of organic matter degradation. 
However, successful application of this model requires knowledge on the 
number of labile groups, their contribution and their reactivity. These are 
unknown and cannot be measured. 

The reactive continuum models do not subjectively partition organic 
matter into a number of pools but consider reactivity as an emergent property 
of the total organic matter pool that continuously declines as reaction pro- 
gresses with time (Middelburg 1989). The two most commonly used models 
are the semi-empirical power model (Middelburg 1989): 


dG 
= = (ONG 3.15 
— = -k(1)G, (3.15) 
where G is the total pool of organic matter and k(t) is time-dependent 
first-order rate constant, and the gamma-distribution reactive continuum 
model (Boudreau and Ruddick 1991): 


ee eee pag 3.16 

F (3.16) 

where G is again the total pool of organic carbon, v is a parameter for the 

shape of the underlying Gamma distribution and km = z% İs the apparent rate 
0 


constant for the mixture: a is a measure for average lifetime of the most 
reactive component and Gp is the initial concentration of organic carbon. 
Moreover, reformulated in a closed form: 


G=G E) (3.17) 


a+t 


it describes the evolution of organic matter as a function of time (Boudreau 
and Ruddick 1991). The power and gamma-type reaction continuum model 
are under certain conditions equivalent and both are particular cases of the 
general q-theory of Bossata and Agren (1995). Although these continuum 
models are less used than the discrete models (Arndt et al. 2013), they have 
been applied for particulate organic matter in soils, sediments and suspended 
particles, phytoplankton degradation experiments and dissolved organic 
matter degradation in lakes. Apparent reactivity constants vary over 8 orders 


51 


52 3 The Return from Organic to Inorganic Carbon 


104 


Photic zone m Sediment traps 


Experiments 102 
« Sediment cores 
— Power Law 


Reactivity (yr) 


10 | Accumulating 104 


Sede es ~~~ ~~~ 


100 
Deep RC CT 10E 
Biosphere a | 104 10? 10° 10? 104 10° 


Observational time window (yr) 


10004 meters below sediment surface 


Fig. 3.9 The reactivity of particulate organic matter is an emergent property that decreases with 
the accumulated time period since it became detritus. Reactivity estimates are based on 
phytoplankton degradation experiments, settling phytoplankton and sediments in the accumulated 
zone. A power law (k = 0.21t~°°85) relates k the first-order kinetic coefficient (yr ') with time 
(t) in yr (Middelburg 1989; Middelburg and Meysman 2007). Fresh phytoplankton initially has a 
high reactivity and this reactivity decreases during degradation because of preferential degradation, 
secondary production and mixing with more refractory components 


of magnitude (Fig. 3.9), and we unfortunately lack analytical protocols to 
measure it directly or indirectly via a proxy. As a consequence, organic matter 
reactivity is a poorly constrained parameter in all biogeochemical models. 


Box 3.2: Carbon consumption in the coastal ocean 

In coastal zone, organic carbon is not only newly produced by phytoplankton, 
but also by microphytobenthos, macroalgae, seagrass, marshes and man- 
groves and imported from terrestrial ecosystems. Duarte and Cebrian (1996) 
presented the fate of autotrophic carbon based on an extensive cross-system 
survey (Fig. 3.10). Autotrophic respiration represents a loss of 26-35% for 
phytoplankton and microphytobenthos, while marine vegetations respire 
between 51 and 69% of their gross primary production. Herbivory losses are 
much higher for phytoplankton and microphytoplankton (26-37%), than for 
seagrass, marshes and mangroves (4—10%) with macroalgae in between 
(~ 16%). Detritus production and respiration were similar among these 
communities (15-27% of GPP) with the balance of organic material being 
exported to adjacent systems. These cross-system patterns in autotrophic 
respiration can be partly attributed to differences in biomass relative to pro- 
duction, because autotrophic respiration scales with biomass (Soetaert and 
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Fig. 3.10 Fate of primary produced materials for marine communities. Blue: autotrophic 
respiration, orange: herbivory; grey: degradation; yellow: export & accumulation (Duarte and 
Cebrian 1996) 


Herman 2009). Differences in herbivory losses between phytoplankton and 
microphytobenthos on the one hand and higher plants on the other can be 
attributed to the general decline in growth rate and palatability and nutrient 
availability from microalgae to higher plants (Duarte and Cebrian 1996). The 
palatability of phytoplankton and microphytobenthos is generally much 
higher than that of terrestrial and littoral plants (mangroves, salt marshes) 
because the latter have lower nutrient and protein contents, higher proportions 
of structural components such as carbohydrates and lignins needed for 
rigidity and sometimes chemical compounds for defences. Grazers on sea- 
grass and salt-marsh plants consequently have different feeding and digestion 
strategies than those on phytoplankton and microphytobenthos. While 
waterfowl, turtles and mammals such as the dugong are iconic grazers of 
seagrasses, many other organisms such as gastropods, isopods and sea 
urchins feed on epiphytes, the microbes growing on the leaves of seagrass. 
These epiphytes function more alike phytoplankton and microphytobenthos 
in terms of productivity, turn-over and their fate. 

Herbivory represents a major loss pathway for phytoplankton and involves 
zooplankton as well as benthic suspension feeders. While benthic and pelagic 
grazers have many commonalities, there are also distinct differences: benthic 
suspension feeders are normally sessile, while planktonic float; benthic sus- 
pension feeders are usually larger and longer-lived with the consequence that 
they overwinter and already have a high biomass before the phytoplankton 
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spring bloom. Moreover, benthic suspension feeders have only access to 
phytoplankton in the lower part of the water column and vertical mixing is an 
important factor governing their growth and their impact on phytoplankton 
dynamics (Herman et al. 1999). 

Coastal systems exchange organisms and detrital material with adjacent 
ecosystems and receive detrital resources via rivers. This additional external 
organic matter is consumed and, as a consequence, the local balance between 
autotrophic and heterotrophic processes is disturbed, with net ecosystem 
heterotrophy as a result (Smith and Hollibaugh 1993). Heterotrophic 
ecosystems often show secondary production, similar to or higher than pri- 
mary production, a dominance of microbial processes, a depletion in dis- 
solved oxygen, and high carbon dioxide levels (Heip et al. 1995). 
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Carbon Processing at the Seafloor 


Most of the global marine primary production (~55 Pg C y_') is consumed in the 
surface, sunlit ocean (~44 Pg C y ') and ocean interior (~8 Pg C y ') and only 
the remaining 2-3 Pg C y ' reaches the sediment floor. Globally most of the 
organic carbon delivered to sediments (~ 90%) is degraded because organic carbon 
burial is low (~0.2-0.4 Pg C y"') and mainly occurs in rapidly accumulating 
coastal sediments. 

In this chapter, we will first discuss organic matter delivery to the sediments, 
then the processes and organisms involved in organic matter degradation and close 
with a discussion of factors governing organic carbon burial. 


4.1 Organic Matter Supply to Sediments 


In Chap. 3 we have seen that organic carbon deposition on the seafloor depends 
primarily on water depth (Fig. 3.4), as a direct consequence of organic matter 
production in the euphotic zone and degradation of detritus during transit from the 
sunlit layer to the bottom of the ocean with more time for degradation in deep-water 
settings than in coastal settings. Consequently, deep-sea sediments receive little 
organic matter, while coastal sediments receive much. This water depth dependency 
of organic matter delivery rates is articulated by lateral differences in ocean primary 
production (production in shallow coastal systems is generally higher than that in 
the open ocean). Marine sediments are therefore often considered donor-controlled, 
i.e. the organic matter consuming benthic communities have no control over their 
food resources. While such a donor-controlled view applies to most deep-sea 
environments, there are other carbon delivery routes for sediments within the 
euphotic zone and for sediments inhabited by animals and these carbon transfers to 
sediments are (partly) mediated by the consuming communities, i.e. under con- 
sumer control (Fig. 4.1; Middelburg 2018). 
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Fig. 4.1 Organic matter supply to sediments (Middelburg 2018). (1) Organic matter settling from 
the water column is deposited at seafloor (donor control; fixed flux upper boundary condition). 
(2) Sediments in the photic zone are inhabited by benthic microalgae that produce new organic 
matter in situ and grazing animals can impact the growth of these primary producers. 
(3) Bioturbating animals transfer labile carbon from the sediment surface layer to deeper layers 
in the sediments. (Vertical axis is depth; horizontal axis is concentration, see Fig. 4.2) 
(4) Suspension-feeding organisms enhance the transfer of suspended particulate matter from the 
water column to the sediments (biodeposition). (5) Sponge consume dissolved organic carbon and 
produce cellular debris that can be consumed by benthic organisms (i.e., the sponge loop) 


Intertidal sediments and coastal sediments within the euphotic zone may support 
microbial photoautotrophs (microphytobenthos). Gattuso et al. (2006) showed that 
globally about 1/3 of sediments within the coastal ocean receive sufficient light for 
photosynthesis. This primary production at the seafloor represent a source 
of ~ 0.32 Pg C yr ' that is not only directly available for consumption, but it is also 
of high quality, because there has been no pre-depositional processing of organic 
matter. This benthic primary production input thus supports about 10-15% of 
global sediment respiration. Moreover, the grazing by micro, meio and macrofauna 
controls the biomass of the primary producers and in this way grazing animals 
control food supply for themselves, as well as for microbes living on the exudates 
and detritus of microphytobenthos (Middelburg et al. 2000; Middelburg 2018). 

Macrophytes such as seagrasses in the subtidal euphotic zone and marsh plants 
and mangroves in the intertidal zone enhance carbon inputs to sediments in two 
ways. They locally produce new organic matter that is used for respiration, invested 
in leaves and stems, allocated to below-ground tissues for nutrient (and water) 
uptake, anchoring and storage. All this tissue can be consumed by microbes and 
animals, although the degradability of structural carbohydrate rich macrophyte 
material is often lower than that of protein-rich microphytobenthos and phyto- 
plankton material (Duarte and Cebrian 1996; Fig. 3.10). Macrophyte communities 
stabilize the sediments and impact the local hydrodynamics with the consequence 
that suspended particles, including particulate organic matter, are trapped within 
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their canopies. This trapped organic matter is often derived from adjacent ecosys- 
tems or phytoplankton from within the system and makes a substantial contribution 
to sediment organic carbon pools in salt marshes (Middelburg et al. 1997) and 
seagrass meadows (Gacia et al. 2002; Kennedy et al. 2010). Moreover, the trapped 
phytodetritus often is more readily degradable than locally produced macrophyte 
material (Barron et al. 2006). 

Suspension feedings worms, bivalves, corals and sponges consume particulate 
organic matter suspended in the water column. Part of the carbon is used for growth 
and maintenance, another part is excreted as faeces and pseudofaeces, which in turn 
are available for microbes and animals in the underlying and nearby sediments. 
Some organisms passively filter water, whereas others actively filter water to obtain 
their resources. Clearance rates (m? of water cleared from food per unit time) can be 
very high, and grazing by benthic suspension feeders often represents the largest 
loss term of phytoplankton in tidal systems (Heip et al. 1995). While some sus- 
pension feeders are spatially distributed rather homogenously, many others (e.g. 
oysters, mussels, coral reefs) occur in specific areas at high densities because of 
local favourable hydrodynamic conditions or positive feedbacks between commu- 
nity density, performance and local hydrodynamics (Herman et al. 1999). These 
suspension feeders locally enhance carbon input to underlying sediments relative to 
that delivered via sinking because they focus laterally available particulate organic 
carbon. Some of the highest carbon inputs to sediments have been reported for 
mussel and oyster beds. Moreover, benthic suspension feeders within the euphotic 
zone directly impact phytoplankton growth and dynamics (Koseff et al. 1993). 

Sandy sediments cover about 50% of the continental shelf seafloor and about 4% 
of the total ocean floor. Their high permeability allows flow of water through the 
connected pore network of these sediments. These flows are driven by current or 
wave induced pressure differences, which may be enhanced through biological 
structures. Active flows of seawater with suspended particles result in additional 
transport of particulate organic matter to sediments. Part of the suspended partic- 
ulate matter entering into these sediments with the water is trapped or utilized by 
organisms (Huettel et al. 2014), thus supplementing the particulate organic matter 
supply to these sediments. This additional organic carbon is highly reactive, and 
particulate organic carbon stocks in permeable sediments are very low (Boudreau 
et al. 2001). 

Another organic carbon provision pathway to sediments involves organisms that 
consume dissolved organic carbon and use it to form particulate organic carbon 
(biomass and excretion products), which can be used by other consumers living at 
the sediment-water interface or within the sediments. This mechanism has been 
documented for encrusting sponges that consume large quantities of dissolved 
organic carbon and rapidly convert this into detritus for benthic consumers (i.e. 
sponge loop, de Goeij et al. 2013). 

It is important to realize that some of these additional pathways involve truly 
new organic carbon (benthic primary production), while the others either convert 
widely available, but poorly useable, dissolved organic carbon into a particulate 
pool (sponge loop) or focus particulate organic matter delivered to a wide area for 
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deposition into small areas inhabited by suspension feeding animals or plants. The 
latter category locally enhances organic matter decomposition at the expense of 
adjacent areas or water-column degradation. Finally, even in donor-controlled 
sediments, animals enhance the transfer of labile organic matter to subsurface 
microbial communities by bioturbation (Box 4.1) and in that way stimulate the 
overall degradation of organic matter (Middelburg 2018). This mechanism will be 
presented in detail below; we will first introduce the consuming organisms. 


4.2 The Consumers 


Organic matter delivered to the sediments represents the major resource for the 
heterotrophic consumers living in the sediments. These organic carbon consumers 
are usually partitioned into size classes: microbenthos (<32 um), meiobenthos (32 
um-0.5/1 mm) and macrobenthos (>0.5/1 mm), the actual size class division 
depending on the environment (smaller for deep-sea than coastal studies). The 
microbenthos includes small eukaryotes (Protists), but is by far dominated by 
prokaryotes (Bacteria and Archaea). The meiobenthos is often dominated by 
nematodes, but foraminifera (protists) also contribute. The macrobenthos comprises 
many different animal groups. The number of organisms declines with size: from 
10° bacteria per ml of sediments (i.e. 10'* cells in the top 10 cm of 1 m’), 10° 
meiofauna m ~ to 10* macrofauna m`? in estuarine sediments (Herman et al. 1999). 
The distances between individuals is thus of the order 10 um, 1 mm and 1 cm for 
microbes, meiofauna and macrofauna, respectively. Microbes are much smaller 
than typical pore sizes of silty and sandy sediments, meiofauna can live within these 
pores, while macrofauna are much larger and consequently move particles when 
active. The mixing of particles and solutes by animals living in sediment is called 
bioturbation (Box 4.1) and has major consequences for carbon processing in sed- 
iments (Meysman et al. 2006; Middelburg 2018). 

The organic carbon consumed by benthic heterotrophs is used for maintenance, 
growth and reproduction. The relative importance of the different size classes to 
carbon consumption scales more or less with their contribution to biomass. In 
coastal sediments, macrofauna typically accounts for 10-25% of the total biomass 
and respiration, with microbes accounting for the remaining, except for a small (few 
%) contribution of protists and meiofauna. Macrofauna number and size (and thus 
biomass) decline stronglywith water depth compared to those of meiofauna and in 
particular microbes (Snelgrove et al. 2018). Carbon consumption in deep-sea 
sediments is almost entirely due to microbes. Accordingly, organic matter con- 
sumption is normally attributed to sedimentary microbes, and the consumption of 
organic matter by animals is therefore ignored in most biogeochemical studies 
(Middelburg 2018). 
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4.3 Organic Carbon Degradation in Sediments 


Organic carbon degradation in sediments is a highly efficient process because most 
of the incoming organic carbon is degraded and only a small fraction eventually 
escapes mineralization to be buried with the accumulating sediments. This fraction 
is often quantified in terms of the burial efficiency (BE): 


_ Fg Fe 


BE = =_—8 _ 
Fo (Fp +R) 


(4.1) 


where Fc is the incoming carbon flux, Fg is the rate of carbon burial and R is the 
total mineralization rate. Burial efficiencies are often presented as a percentage and 
range from a fraction of a percent in deep-sea sediments up to tens of a percent in 
rapidly accumulating coastal sediments (Canfield 1994; Aller 2013). In most 
oceanic sediments, R >> Fp, and total respiration is a reliable proxy for quantifying 
the organic carbon delivery rate to sediments (Fc). 

This high efficiency of organic matter degradation within in sediments might 
seem paradoxical because sediments receive organic matter that has already been 
processed extensively before deposition. Pre-depositional processing of organic 
matter not only lowers the amount of organic matter delivered, but also the quality 
because of the preferential use of labile organic carbon by organisms (Fig. 3.9). 
However, when particulate organic matter is transferred from the water column, 
with a typical residence time of weeks to months, to the top layer of sediments, the 
time available for degradation changes by orders of magnitude (10°-104 years). 
This increase in time available for processing organic matter more than compen- 
sates for the loss of reactivity during particle settling; hence, organic matter 
degradation resumes. 

To understand the dynamics and distribution of organic matter in sediments we 
will again use the simple diffusion-advection-reaction model introduced in Box 1.1. 
Modelling the distribution and dynamics of organic carbon in sediments is rather 
challenging because the transport of organic carbon, as well as its production and 
consumption, are biologically controlled, and we cannot resort to physics-based 
laws for transport. Moreover, we are dealing with a two-phase problem (particles 
and water). The first simplification is that we only consider particulate organic 
carbon. For a sediment with a porosity (volume water/total volume) of 0.7, a 
particulate organic carbon content of 1 wt% and a very high pore-water dissolved 
organic carbon concentration (1 mM), dissolved organic carbon contributes about 
1%o to the sediment carbon stock. Moreover, the reactivity of dissolved organic 
carbon is relative low. Two, gradients in porosity due to compaction and biotur- 
bation are ignored. Three, organic carbon degradation follows first-order kinetics 
(one-G model, Berner 1964; Box 3.1). Four, electron acceptors are abundant and 
their nature does not impact organic carbon degradation. Five, solid-phase organic 
carbon is transported by bioturbation (Box 4.1) and by the net accretion of 
sediment. 
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Transport of particulate organic carbon due to animal activity is the dominant 
transport process in sediment inhabited by animals. The transport resulting from 
bioturbation can usually be well described by diffusion (Goldberg and Koide 1962; 
Berner 1980), in particular when the benthic community is diverse (Meysman et al. 
2006, 2010). It is only in sediments where chemical conditions (e.g. lack of oxygen, 
pollution) or physical conditions (e.g. unstable sediments that are repetitively dis- 
turbed by waves/tides) restrict macrofauna that bioturbation can be ignored and 
organic carbon transport due to sediment accretion dominates. Finally by assuming 
steady-state, the resulting equation for organic carbon (G) is: 


LG dG 
=D ye 4.2 
ba ds a 


0 
where G is labile particulate organic carbon, x is depth (positive downwards), Dy is 
the bioturbation coefficient (cm? yr’), w is the sediment accumulation rate (cm 
y ') and k is a first-order rate constant (yr '). Considering that all labile organic 
carbon is eventually consumed (at depth), we can solve the equation for specific 
conditions at the sediment-water interface (x = 0). It is instructive to distinguish 
between sediments that are strictly donor-controlled (i.e. a fixed flux upper 
boundary) and sediments with a fixed concentration in the top layer. The former 
applies to most sediments, while the latter would represent coastal sediments with 
multiple deposition—resuspension events and lateral particulate organic transport 
pathways (Rice and Rhoads 1989). If we know the concentration of organic carbon 
at the sediment-water interface (Go), the solution is (see Box 1.1): 


— Vw? +4kD 
G= Ge" with a= eov (4.3) 
b 


Alternatively, if the flux of organic carbon delivered to the sediment (F) is known 
the solution is: 


G= - mm (4.4) 
= —Dia-ew f ` 


where o is the same as in Eq. 4.3. 

Figure 4.2 shows the distribution of organic carbon as a function of sediment 
depth, for a fixed reactivity and sediment accumulation rate, and for different 
bioturbation coefficients. The distribution of organic carbon declines exponentially 


k 


with depth and the attenuation is to a first-approximation governed by 4/5; i.e. 


more mixing (higher D,) will flatten organic carbon concentration versus depth 
profiles (because the contribution of w to the attenuation coefficient « is very small). 
In other words, particle mixing by moving animals will transfer organic carbon to 
larger depths. This is articulated in sediments with a fixed concentration in the top 
layer. For these sediments bioturbation will increase the total amount of organic 
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Fig. 4.2 Organic carbon concentration versus depth profiles for different bioturbation coefficients 
(0.1, 1, 10, 100 cm? yr) for a fixed concentration of 20 mg C gr! (left) and a fixed flux of 
36 mg Ccm * yr | (right). The imposed flux has been tuned to generate an identical concentration 


at the sediment-water interface for a Dp of 10cm? yr. Model parameters: k = 1 yr‘; 


w=0.1 cm yr! 


matter available for the benthic community. In other words, animals have, to a 
certain extent, control over their food supply. The other major factor governing 
organic matter profiles in sediments is the reactivity of organic matter. More 
reactive organic matter is degraded shallower and does not penetrate as deep into 
the sediment as refractory organic matter (Fig. 4.3). Moreover, the inventory of 
organic carbon in surface sediments is higher when the reactivity of organic matter 
is lower, in particular for a fixed amount of carbon supplied per unit time. 
Evidently, when applying this type of models, it is important to select the 
appropriate upper boundary conditions. While deep-sea sediments can best be 
described by an organic carbon flux boundary condition, the proper upper boundary 
condition for coastal sediments is not clear. A fixed concentration upper boundary 
conditions may be more appropriate for shelf systems in which lateral organic 
carbon fluxes dominate, e.g. due to a combination of resuspension-deposition 
cycles and residual tidal flow. Animals collecting organic matter in the surface 
sediment layer and bioturbating it to deeper layers make thus additional organic 
matter of high quality available to microbes. Moreover, first-order reactivity 
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Fig. 4.3 Organic carbon concentration versus depth profiles for different first-order reactivity 
coefficients (0.1, 1, 10, 100 yr‘) for a fixed concentration of 20 mg C gr! (left) and a fixed flux 
of 36 mg C cm”? yr | (right). The imposed flux has been tuned to generate an identical 
concentration at the sediment-water interface for k = 1 yr '. Model parameters: Dp = 10 cm? 
yr! 

constants (k) vary over orders of magnitude (Fig. 3.9) and unfortunately cannot be 
directly measured; they have to derived by fitting a model to observations. This 
approach is problematic for a number of reasons. One, fitting of a model to 
observed organic carbon versus depth profiles allows constraining the depth 


E but not k or D» individually (Soetaert et al. 1998; an 


;w=0.1 cm yr" 


attenuation coefficient 
iconic example of an identifiability problem). Two, both the organic matter reac- 
tivity and bioturbation coefficient may depend on depth within sediment. The 
former because organic matter at depth has on average been exposed longer to 
degradation and the latter due to the depth distribution of infauna that results in a 
depth gradient in mixing intensity. 

The depth of the bioturbated zone (L), the reactivity of organic carbon (k) and the 
mixing intensity (D,) are all related (Boudreau 1998). If we assume that the depth of 
bioturbation is governed by the availability of food for animals and set a lower limit 
for organic carbon at depth L (G) at 1% of the organic carbon at the surface (Go): 
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G = 
a = 0.01 = eV at (4.5) 


0 


Making use of In(0.01) ~ 4.6, this equation can be re-arranged to 


L= 4.6/7 (4.6) 


The depth of the mixed layer is thus directly related the reactivity of the organic 
matter and the mixing intensity. This equation is highly similar to the more elab- 


orate result of Boudreau (1998) L = 4 22 that is based on a resource-feedback 


model. Figure 4.4 shows the depth of the bioturbated zone as a function of organic 
matter reactivity and bioturbation intensity. The observed range of Db and k in 
marine sediments is also indicated. Bioturbation depth in marine sediments vary 
between 1 and 20 cm, with a global average of about 10 cm (Boudreau 1998). This 
can be obtained either by high mixing, if organic matter is highly reactive as in 
coastal sediments, or by low mixing, with less reactive matter in the deep sea. 


4.4 Consequences for Sediment Biogeochemistry 


Degradation of organic carbon to inorganic carbon requires an electron acceptor, e.g. 
oxygen. In the water column, oxygen supply is usually sufficient to accept all electron 
released during respiration, but in sediments solutes such as oxygen are transported 
primarily by molecular diffusion, and oxygen supply rates often do not match oxygen 
consumption rates with the consequence that most sediments become anoxic at some 
depth. In the absence of oxygen, organic matter is either respired anaerobically or 
fermented. Fermentation refers to an internal redistribution of electrons among carbon 
atoms with the result that carbon dioxide and methane are formed. 


2CH20 = CO2 + CH4 


Anaerobic respiration and fermentation of particulate organic matter involves 
multiple steps in a complex network: hydrolysis and fermentation of larger mole- 
cules to smaller molecules and finally the use of alternative electron acceptors 
(Jørgensen 2006). These intermediate products include short-chain fatty acids such 
as formate, acetate, propionate and butyrate, and hydrogen. The latter has a high 
turn-over (minutes) because of rapid production and consumption, e.g. during 
hydrogenotrophic processes such as carbon dioxide reduction to methane (e.g. 
Beulig et al. 2018). 
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Fig. 4.4 Relationship between organic carbon reactivity (k) and particle mixing (Db) for a 
number of bioturbation layer depths (0.1, 1, 2, 5, 10 and 20 cm). Sediments receiving reactive 
organic material are more intensively mixed than those receiving refractory organic matter (blue 
background box). The global average bioturbation depth is about 10 cm (Boudreau 1994) 


The most important alternative terminal electron acceptors are nitrate and nitrite, 
manganese and iron oxides and hydroxides, and sulphate, and they are used in this 
sequence largely because of multiple factors, including energy yield. Table 4.1 
provides the overall reaction for organic carbon respiration based on nitrate, 
manganese oxide, iron oxide and sulphate and the corresponding change in Gibbs 
free energy. 

The Gibbs free energy change declines from aerobic respiration to denitrification 
(nitrate reduction to dinitrogen gas) and metal oxide reduction to sulfate reduction 
because of the formation of reduced compounds (metabolites) such as ammonium, 
manganese (II), iron (II), hydrogen sulphide and methane (Fig. 4.5). These reduced 
metabolites contain a substantial amount of energy that was originally part of the 
organic matter. Once produced these reduced metabolites can react (after diffusion 
upward) with an electron acceptor (higher in the redox cascade). This re-oxidation 
of reduced metabolites guarantees that sediments efficiently utilize all organic 
matter energy delivered (Fig. 4.5). For instance, methane generated at depth can be 
oxidized anaerobically with sulphate, metal oxides and nitrate/nitrite and aerobi- 
cally. Some of the metabolites generated (e.g. sulphide) may react with sedimentary 
iron containing minerals. The iron sulphide formed may be buried, but the majority 
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Table 4.1 Standard gibbs free energy change of organic matter degradation pathways (based on 
Berner 1980) 


Degradation pathway Reaction AG° 
(kJ mot! | 
of CH,0) 
Aerobic respiration CH,0 + O; = CO, + H2O —475 
Denitrification CH,0 + 0.8 NO; <= 0.4 N2 + 0.8 HCO; + 0.2 —448 
CO2 + 0.6 H2O 
Manganese oxide CH,O + 3 CO, + 2 MnO, + H2O & 2 Mn” +4 —349 
reduction HCO; 
Iron(hydr)oxide CH,0 + 7 CO, + 4 Fe(OH); = 4 Fe™ +8 —114 
reduction HCO; +3 H,O 
Sulfate reduction CH20 + 0.5 S04” © 0.5 H2S + HCO; =I]. 
Methanogenesis CH20 = 0.5 CO; + 0.5 CH4 —58 


The decrease in Gibbs free energy of the reaction results from the accumulation of energy in the 
reaction products, i.e. the metabolites 


is usually re-oxidized after bioturbation to surface layers where oxygen or nitrate 
are present (Jørgensen 1977; Berner and Westrich 1985). Oxidation of these 
reduced metabolites often involves multiple reactions: e.g. hydrogen sulphide is 
oxidized to intermediate compounds such as elemental sulfur, thiosulfate and sul- 
phite, which may then be oxidized further or disproportionated (i.e. splitted into a 
more reduced and more oxidized compound). 

This efficient re-oxidation of reduced metabolites has two important conse- 
quences for studying the sedimentary carbon cycle. First, the energy released during 
re-oxidation reactions can be used to fix inorganic carbon, i.e. it supports 
chemolithoautotrophic processes in sediments. Chemolithoautotrophy in marine 
sediments is still poorly constrained, but likely on the order of 0.3-0.4 Gt C yr! 
(Middelburg 2011). Secondly, total oxygen consumption is not only linked to 
aerobic respiration, but also includes oxygen consumption during re-oxidation 
processes. In organic-carbon limited sediments in the deep-sea, oxygen supply by 
diffusion is sufficient and all organic matter mineralization occurs aerobically, and 
oxygen use is primarily due to oxidation of organic carbon to carbon dioxide and 
nitrification (oxidation of ammonium to nitrate). In sediment with high carbon 
loadings, most of the oxygen consumption is due to the oxidation of reduced sulfur, 
iron, manganese and methane, and total oxygen consumption provides an excellent 
proxy for total sediment respiration irrespective of actual respiration occurred 
aerobically or anaerobically (Fig. 4.6). In fact, total respiration of sediments is 
equal to total oxygen consumption minus the burial of reduced iron and sulphide 
and the escape of nitrogen gas over the sediment-water interface. These contribu- 
tions are usually limited to a few percent and oxygen consumption can, therefore, 
be used as a measurable proxy for total respiration without the need to study all 
microbiological reduction and oxidation processes in detail. 

Oxygen distributions in sediments are normally modelled with a simple model 
balancing the consumption of oxygen with the diffusive supply: 
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Fig. 4.5 Conceptual model of organic matter (OM) degradation and reoxidation pathways (based 
on Jørgensen 2006, and Middelburg and Levin 2009). The red arrows reflect the fate of 
(energy-rich) substrates released during anaerobic mineralization 
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Fig. 4.6 Conceptual model of tight coupling sulfate utilization for respiration of organic carbon 
and sulphide oxidation by oxygen, with the net result that one mole of oxygen use corresponds 
with one mole of organic carbon oxidized to carbon dioxide (after Aller 1994) 
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dO 

D; Pe = Ro (4.7) 
where Op is the oxygen concentration, D, is the diffusion coefficient of oxygen in 
sediments and Rog is the consumption of oxygen. Oxygen consumption rates can be 
described with a depth and oxygen concentration independent rate because 
re-oxidation is the dominant oxygen sink. The solution for this zero-order kinetics 
equation with a fixed concentration at the sediment-water interface (Co) is 
(Bouldin 1968): 


R 2CoR x\2 
Ca" %2 ( re) + C= a17) (4.8) 


where 6 is the oxygen penetration depth: 


_ PCD, 


This relation shows that the oxygen penetration depth is inversely related to the total 
respiration. High respiration rates in coastal sediment result in shallow oxygen pen- 
etration (mm scale), while low respiration rates in organic carbon starved deep-sea 
sediments allow oxygen penetration to meters (Fig. 4.7). Moreover, oxygen pene- 
tration depth exhibits a square-root dependence on bottom-water oxygen levels (Co). 
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Fig. 4.7 The relationship between oxygen consumption and oxygen penetration depth. High 
fluxes of oxygen into sediments imply shallow penetration. Oxygen penetration varies 
from <1 mm in coastal sediments to >1 m in deep-sea sediment receiving low carbon loadings 
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4.5 Factors Governing Organic Carbon Burial 


Following active carbon degradation in the water column and continued mineral- 
ization in surface sediments, at some depth in the sediment all labile organic matter 
has been utilized and a background level of organic carbon remains (Figs. 4.2 and 
4.3). This organic carbon is buried down with the accretion of sediments. Organic 
carbon burial is, thus, the product of sediment accumulation and organic carbon 
content: 


Fs = w(1 — 6)p.G (4.10) 


where w is the sediment accumulation rate (cm yr '), ¢ is the porosity, p, is the dry 
density (gr cm *), G is the concentration of organic carbon (g C gr '), with Fg 
expressed in g C cm * yr". It is important to realize that w refers to the long-term 
accumulation rate at depth and not the rate of sediment deposition at the surface. 
The concentration of organic carbon should theoretically be at the depth where no 
further degradation occurs. Such a depth would not exist if organic matter decay 
follows a reactive continuum: organic matter is reactive on all timescale. However, 
in most sediments organic carbon concentration rapidly declines with depth and 
then remains almost invariant. The precise reference depth is not so important as 
long the gradient can be neglected. 

Differences in organic carbon burial are primarily determined by differences in 
sediment accumulation and carbon contents, and only secondarily by differences in 
porosity and density. Sediment porosities vary 0.4—0.8 and dry sediment densities 
from 2.6 gr cm ° for mineral to ~2.0 gr em ° for peaty sediment. Organic carbon 
concentrations depend on the carbon input to sediments and the extent of decay 
before burial. Any process that enhances organic carbon input or impedes degra- 
dation elevates organic carbon concentration at depth and consequently carbon 
burial. These two factors relate directly to the old discussion on whether high export 
production (carbon delivery to sediment) or an oxygen related effect on preservation 
governs burial of organic carbon. The idea is that organic matter degradation in the 
presence of oxygen is more efficient than without (Box 4.2). 

The by far most important factor governing organic carbon burial is the delivery 
of inorganic particles that results in sediment accretion. Sediment accumulation in 
the ocean vary about four to five orders of magnitude, from less than 0.1 cm per 
thousand year in the deep-sea to centimetres per year in rapidly accumulating 
coastal sediments (Aller 2013). Most organic carbon burial, therefore, occurs in 
coastal sediments. An increase in inorganic particle deposition should, in principle, 
not increase carbon burial, because the organic carbon will just be diluted, but 
observations show otherwise. In nature, particulate organic carbon is mainly 
associated with mineral surfaces and the higher the accumulation of minerals, the 
more mineral surfaces, with associated organic matter, are buried. Vegetated coastal 
sediments are another site of high organic carbon burial for two reasons. First, the 
presence of macrophytes stabilizes sediment and enhances deposition via canopy 
trapping of particles. Secondly, the net ecosystem carbon balance of these 
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communities results in the accumulation of plant derived organic matter. The fol- 
lowing flux balance can be used to investigate the impact of sediment accretion and 
net ecosystem production on the burial of organic carbon in vegetated sediments: 


w(1 _ P) PsCexternal ar PorganicCtoM NEP = w(1 T P)PsCburied (4.11) 


where w is sediment accumulation rate (cm yr '), ọ is the porosity, p, is the dry 
density of bulk sediment (gr cm °), Cexternal is the concentration of organic carbon 
(gC gr ') in deposited sediments remaining after degradation, Porganic 18 the density 
of plant organic material, CtoM converts organic C to organic matter, and Cpuriea iS 
the concentration of buried organic matter. Primary production in cordgrass habitats 
is very high (~2 kg C m° yr ‘) with net carbon accretion of ~0.3 kg Cm? yr ! 
(Middelburg et al. 1997). If all that cordgrass detritus is retained in the system, 
marsh sediment accretion would be ~0.5 mm per year, not enough for keeping 
pace with accelerated global sea-level rise (>2 mm per year) and a mineral starved, 
peaty marsh may therefore drown. The combined effect of enhanced particle 
trapping and retention and macrophyte derived organic matter are the reason that 
vegetated coastal sediments represent a major term in global marine carbon burial, 
i.e. blue carbon (Duarte 2017). 


Box 4.1: Bioturbation 

Bioturbation refers to the reworking of sediments (particles and water) due 
the activities of organisms, in particular animals (Meysman et al. 2006). This 
reworking has many consequences, including the creation of a 
three-dimensional, heterogeneous distribution of sediment properties. In fact, 
this animal induced heterogeneity is clear from sediments deposited during 
the last 0.5 billion years and palaeontologists have therefore named this 
period Phanerozoic, after the old greek words phanerds (visible) and zōé 
(life). Moreover, the small-scale heterogeneity stimulates biogeochemical and 
biological diversity. It also has consequences for porosity, permeability, 
sediment stability and roughness. 

Besides these effects on sediment properties and characteristics, animal 
activities have a major impact on the transport of particles and solutes. Par- 
ticle mixing, bioturbation sensu stricto, occurs through the construction and 
maintenance of structures such as burrows and the movement of particles 
while feeding (ingestion and defecation). Particles are moved vertically and 
horizontally, but the consequences are usually more prominent vertically 
because most sediment characteristics show strong vertical gradients. The 
particles dispersed may be reactive or unreactive minerals, organic matter as 
well as organisms and their cysts, eggs and remains. Many animals construct 
and live in burrows and ventilate these with overlying oxygenated water for 
respiration or food supply. In permeable sediments, this burrow water can 
enter the sediments and in permeable and non-permeable sediments diffusive 
exchange occurs between burrows and adjacent sediments because of 
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concentration gradients. Many animals are involved in particle mixing and 
ventilation. Particle mixing is pivotal for optimal functioning of sediments: 
organic matter and solid-phase oxidants such as iron oxides are mixed 
downwards, while reduced substance such as iron sulphides are mixed 
upwards. Bio-irrigation enhances exchange of oxygen, nutrients and other 
substances between water and sediment. 

Particle mixing is the result of multiple activities of a diverse assemblage 
of large animals. While the behaviour of certain species may cause a distinct 
particle transport pattern, the net result of many particle mixing events can 
normally be described very well as diffusion. Bioturbation is then quantified 
in terms of a biodiffusion coefficient for particles (D,; cm? yr’). This 
biodiffusion coefficient can be derived from modelling concentration versus 
depth profiles of a tracer A with a known reactivity (e.g. radioactive decay). 
Equation 4.12 balances transport due to sediment accumulation and biodif- 
fusion with radioactive decay at steady state: 


PA dA 
w kA (4.12) 


aes a dx? dx 


where w is the sediment accumulation rate (cm yr’) and k is the radioactive 
decay constant (yr '). These macroscopic empirical D, values can be 
decomposed into discrete particle properties step length (à) and the period 
between two displacements (t) (Meysman et al. 2010): 


Pa (4.13) 


The step length of particle displacement scales with organism size, and it is 
for this reason that meiofauna and similarly sized organisms contribute little 
to particle displacement. 


Box 4.2: Oxygen and Organic Matter Preservation 

The impact of oxygen on organic matter degradation and preservation has 
been studied extensively because of its relevance for the formation of oil and 
gas source rocks, the reconstruction of paleoenvironments, and the projection 
of carbon cycling in the future warmer, low oxygen ocean. Although most of 
the evidence is correlative, there is now consensus that more organic carbon 
is preserved under low oxygen conditions (Middelburg and Levin 2009). 
Experimental studies have shown that oxygen has little if any impact on 
microbial organic matter degradation, but is needed for the degradation of the 
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most refractory fraction of sediment organic carbon (Hulthe et al. 1998). 
Natural experiments provided by turbidite deposition and cyclic 
de-oxygenation events have revealed that organic matter preserved under 
anoxic conditions can be re-activated when exposed to dissolved oxygen 
(Moodley et al. 2005). Bottom-water oxygen levels have been shown to 
increase the carbon burial efficiency and the amount of carbon preserved per 
unit area of reactive surface area (Canfield 1994; Hartnett et al. 1998). Lack 
of oxygen not only increases the quantity of organic carbon buried, but also 
its nature. Organic carbon buried under anoxic conditions is usually less 
degraded and that is reflected in organic matter proxies (see Chap. 6), such as 
hydrogen index and the amino-acid degradation index (Middelburg and Levin 
2009). Elevated organic carbon burial under anoxic conditions is a negative 
feedback in Earth System dynamics. 

Although the impact of dissolved oxygen on carbon preservation has been 
well documented, and we have made progress studying the implications, 
there is still little understanding of why there is more organic carbon burial 
under anoxic conditions. Most organic matter is buried in ocean margin and 
coastal sediments that are anoxic below a few mm to cm, irrespective of the 
presence of oxygen in the bottom waters. This implies that the dissolved 
oxygen effect should be either pre-depositional or related to the changes in the 
benthic community processing the delivered organic matter. There is evi- 
dence for both. Organic carbon flux attenuation in an anoxic water column is 
less than that in an oxic water column, perhaps due to the absence of zoo- 
plankton under anoxic conditions, and the organic matter delivered to anoxic 
sediments appears to be less reactive towards degradation (Keil et al. 2016). 
Changes in bottom-water oxygen levels have consequences for the benthic 
community composition because of the decrease or disappearance of animals 
(Jessen et al. 2017). Benthic animals play a major role in sediment ecosystem 
functioning, through their interactions with microbes: particle mixing delivers 
labile organic carbon to microbes living at depth and bio-irrigation provides 
microbes with fluctuating oxygen levels, so that they can efficiently process 
organic matter (Aller 2013). 
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Biogeochemical Processes 5 
and Inorganic Carbon Dynamics 


Organic matter production and degradation result in the consumption and release of 
carbon dioxide; oxidation and reduction reactions involve the production or con- 
sumption of protons; carbonate mineral formation and dissolution reactions cause 
consumption or release of carbonate and bicarbonate. All these processes have the 
potential to increase or decrease proton concentrations, i.e. change the pH, yet 
marine pH varies over a rather small range (about one unit only). This is due to the 
buffering of seawater. For a complete understanding of the impact of carbon dioxide 
on organisms, the role of organisms in carbon dioxide dynamics, and the role of 
marine systems in climate change (e.g., carbon dioxide uptake and ocean acidifi- 
cation; Box 5.1), we have to understand the chemistry of the carbon dioxide system 
in water and for this we need to understand what governs pH dynamics. 

In this chapter, we will first refresh the basics of inorganic carbon chemistry in 
water and seawater, then introduce alkalinity and buffering, and discuss how bio- 
logical processes and mineral formation are impacted and impact carbon dioxide in 
marine systems. 


5.1 The Basics 
Pure water can dissociate into protons and hydroxide ions: 
H20 + H* + OH” (5.1) 


This reaction occurs virtually immediately and one can thus assume equilibrium 
between the three species ((H2O], [H*], [OH ]): 


[H ][0H ] 


K= moj 


or K, = [H+ ][OH7] (5.2) 
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where K,, denotes the equilibrium constant for water self-ionisation. This type of 
equilibrium relation is also known as a mass-action law. Having two unknowns 
({H*], [OH ]) that are related via the above equilibrium relation, we have one 
degree of freedom: i.e. if we choose [H*], then [OH ] is set to ay 

Next, we consider pure water in equilibrium with a gas containing CO). Part of the 
CO, will remain in the atmosphere (CO,,,)), but part will dissolve in water (COp,aq)). 
This dissolved CO aq) Will react with water to form carbonic acid (H2CO3) according to: 


HO + COraq) > HCO; (5.3) 


This equilibrium is somewhat slow and the equilibrium constant of hydration K, = 


mori is rather small (~ 0.002); almost all dissolved carbon dioxide remains 
in the form of COzaq. Moreover, it is analytically not possible to distinguish 
between COd(aq) and H2CO3(aq); they are therefore lumped together and usually 
termed CO} or HCO}. However, carbonic acid H2CO3,aq) is a moderately weak 
acid (K ~ 10-*°), while the combined HCO} is a weak acid (K ~ 107°). For 
reasons of notational simplicity, we will term it H,CO 3 from now onwards. 

Carbonic acid is a weak diprotic acid and partly dissociates first into a bicar- 
bonate ion (HCO; ) and a proton, and subsequently, the bicarbonate is dissociated 
partly into carbonate ion (CO;” ) and a proton. The relevant reactions are: 


HCO; + HCO; + Ht (5.4) 
HCO; = CO} + Ht (5.5) 


for which we can write equilibrium relations: 


_ [HCo;][H*] 
Wacol = 
and 
_ (cory 
O= Heos] S7 


where K, and K, are the first and second equilibrium constants (107635 and 10 '°% 
in freshwater at 25 °C). The relative concentrations of [H.CO3], [HCO3 ] and 
[CO3* ] are governed by the pH (—logio[H*]) of the solution as depicted by a 
Bjerrum plot (Fig. 5.1). Carbonic acid is the dominant species at pH values below 
the pK, (—log,oK,), bicarbonate dominates between the pK, and pK, values and 
the carbonate ion dominates at pH values above the pK>. 

The carbonic acid in solution and carbon dioxide in the atmosphere are related 
via Henry’s law: 
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Fig. 5.1 Bjerrum plot showing the distribution of carbonic acid (red), bicarbonate (black) and 
carbonate (blue) ions as a function of pH in freshwater (dashed line) and seawater (solid lines). 
DIC = 2000 uM 


[H2CO3] = Ky x pCO, (5.8) 


where Ky is the Henry constant for CO% solubility in water (mol kg! atm™') and 
pCO? is the partial pressure of CO, (atm). The equilibrium constants Ky, Ki, K2 
and Ky, are functions of temperature, pressure and salinity in seawater. 

Accordingly, for the CO,—H20 system, we have five unknown concentrations 
({H»CO3], [HCO; ], [CO;” ], [OH ] and [H*]) and four equilibrium relations: 
water self-ionisation (5.2), Henry’s law (5.8) and the first and second equilibria of 
carbonic acid dissociation (5.6, 5.7). To solve the system, we need one additional 
relation. Natural waters are uncharged, and we can thus use the charge balance 
equation for this system: the positive charge of protons should be balanced by the 
negative charge of hydroxide, bicarbonate and carbonate ions. 


[H+] = [OH] + [HCO;] +2[COz] (5.9) 
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Note that the carbonate ion is counted twice in a charge balance because of its 
double charge. Alternatively, one can define a proton balance equation, a mass 
balance for protons (Butler 1982): 


H7] = H" lno + [H* Jaco, (5.10a) 
or its equivalent 
[H+] = [087] + [HCO;] +2[C0}], (5.10b) 


This proton conservation equation balances excess protons on the left-hand side 
with the recipe on the right-hand side. Proton mass balances are always relative to a 
proton reference level, e.g. Equation (5.10a, b) is relative to HCO; and in this case, 
it is identical to the charge balance. Mathematically, the system is now fully defined 
with 5 unknown species linked via 5 equations. Moreover, the system remains fully 
determined if another component is added for which the total concentration and 
equilibrium distribution among species is known. 


5.2 The Thermodynamic Basis 


So far, we have ignored non-ideal behaviour of gases, liquids and solutes that have 
an impact on equilibria. Although the name equilibrium constant suggests that it is a 
constant, it is in fact not, because of its dependence on temperature, pressure and 
the composition of the solution. This dependence relates to the effect of tempera- 
ture, pressure and mixture composition on the Gibbs Free energy. Formally, at 
equilibrium, the change in Gibbs free energy AG“ (J mol‘) is related to the 
thermodynamic equilibrium constant (Ki"°"™) via AG’ = —RT In Kiherm where T 
is in Kelvin and R is the universal gas constant (8.314 J mol! K`’). 

There are two generic pathways towards obtaining equilibrium constants for 
complex natural solutions. The first approach involves calculating thermodynamic 
constants from the well-documented and tabulated standard Gibbs free energies of 
formation and correcting these for temperature and pressure of the system. The 
thermodynamic equilibrium constant is then expressed in activities of the molecules 
involved, rather than concentrations. For example, equilibrium relation (5.3) for the 
thermodynamic KP®™ should formally be written as: 


{H2CO3} 
where the {} indicate that activities, rather than concentrations [] are used. Con- 
centrations and activities are linked via activity coefficients (y), which account for 
both electrostatic interactions, as well as formation of ion-pairs among the various 
ions in solution: 
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_ {HCO;}{H*} _ [HCO;][H*] „ HCO; H* 
{H2C03} [H2CO3] VH CO3 


therm 
K 1 


These activity coefficients are usually calculated using the Pitzer model (Millero 
2007) that can be applied to highly complex media, such as brines. 

The alternative way of estimating equilibrium constants, used by most chemical 
oceanographers and marine biogeochemists, is to experimentally determine stoi- 
chiometric constants as a function of temperature, pressure and salinity and rep- 
resent this dependency by a polynomial function. The underlying idea is that the 
relative composition of seawater is rather constant, that the reaction as such does 
not change the composition of seawater and that direct measurements in seawater 
provide more accurate results than using thermodynamic data with activity coeffi- 
cient corrections: 


. ___ [ACOs ]|H*| 
Ki rsp) = mco) 


where the * superscripts indicates that it is a stoichiometric quantity, rather than 
thermodynamic constant and the subscripts T, S and P stand for temperature, 
salinity and pressure. The temperature, salinity and pressure dependence of stoi- 
chiometric constants result in rather complex and cumbersome expressions 
(Table 5.1), but there are many computer programs available to facilitate their 
determination. 

Figure 5.2 shows the stoichiometric constant expressed as pK,* as a function of 
temperature and salinity. pK,* values decrease with increasing temperature and 
salinity, or in other words, carbonic acid is more dissociated with increasing tem- 
perature or salinity. The ratio of the stoichiometric equilibrium constant (K,*) at 
salinity 35 to the thermodynamic equilibrium constant (KY) in pure water is 
about 4.4-5. This salinity effect has consequences for the speciation of dissolved 
inorganic carbon in seawater (Fig. 5.1). Similar dependencies apply to other stoi- 
chiometric constants. Moreover, this stoichiometric approach can also be used to 
accurately quantify the solubility of minerals in seawater and will be used from now 
on (while dropping the * superscript). 

The solubility of gases, including COs, is also a function of temperature, salinity 
and pressure (increasing salinity and temperature lower the solubility of gases). 
Thermodynamically, the fugacity of carbon dioxide is linked via the Henry constant 
to dissolved carbon dioxide: [H2 C03] = Ky x fCO2. The fugacity is not exactly the 
same as the partial pressure, the product of CO, mol fraction Gco) times total gas 
pressure, as presented above (in Eq. 5.8). Because most readers are more familiar 
with the pCO», than fCO2, we will use the former notation, noting that they are 
different. 


82 5 Biogeochemical Processes and Inorganic Carbon Dynamics 


Table 5.1 Acid-base reaction in seawater and examples of expression of their stoichiometric 
constants (Dickson et al. 2007). T is absolute temperature (K), S is salinity 


Reaction Equilibrium relations Stoichiometric constant expression 


H0 & H* + OH K, = [H*|[OH-] _ —13847.26 
k T 


In(Ky) + 148.9652 — 23.6521 In(T) 


118.67 
+ i = 5.977 + 1.0495in(7)) S°5 — 0.016155 


co HCO. H: = Ky x pC 100 T 
A Saree FAC) tsa see) In(Ky) =93.4517( — ) — 60.2409 + 23.3585 In( — 
T 100 
if Tae 
+s(o03s0 — 0.023656 cn + 0.0047036 ls) ) 
EN —3633.86 
AEO ECO E | e Eea In(K,) == + 61.2172 — 9.67770 (T) 
+ 0.0115558 — 0.0001 1525? 
HCO; = CO} + Ht e 11778 
3 3 o= ET In(Ka) = — 25.9290 + 3.16967 In(T) 
+ 0.01781 — 0.0001 1228? 
B(OH), + HO = g, — Bonde] —8966.90 — 2890.53595 — 77.9425 + 1.728815 — 0.09965? 
B(OH); + H* B = [B(OH] In(Kg) = 7 


+ (148.0248 + 137.19425°° + 1.621425) 
+ (—24.4344 — 25.0858°° — 0.24748) In(T) + 0.053105S°5T 


5.3 Analytical Parameters of the CO, System 


Not all individual species of the CO, system (pCO2, [H2CO3], [HCO3 ], [CO;” ], 
[OH ] and [H*]) can be measured directly in seawater, and there is also no need to 
because they are interlinked via equilibrium, mass balance and charge conservation 
equations (Dickson 2011). There are five parameters that can be measured: 


(A) The total concentration of dissolved inorganic carbon, often abbreviated as 
DIC, XCO, or Cy: DIC = [H2CO3] + [HCO; ] + [CO ], is normally 
measured by acidifying the sample, stripping the evolved gas and measuring 
the total CO2 content. 

(B) The partial pressure (fugacity) of carbon dioxide (pCO ) can be obtained by 
measuring the gas phase composition in equilibrium with the water and the use 
of Eq. (5.8: Henry’s law). 

(C) The pH can be measured directly using electrodes and by colorimetry using an 
indicator dye (Dickson 2011), which then provides data for the hydrogen 
concentration. However, there is also a large body of research based on the 
NIST scale (National Institute of Standards and Technology, formerly known 
as NBS scale), the total (pH) and the seawater (pHsws) scales. The latter two 
incorporate association of free protons with bisulphate, or bisulphate and 
fluoride, respectively, and are linked to the free scale as follows: pH =~ 
pHr + 0.11 and pH © pHsws + 0.12. 

(D) The carbonate ion concentration (COF via measurement of the ultraviolet 
absorbances of lead carbonate complexes in seawater. 
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Fig. 5.2 The dependence of pK, on the temperature and salinity of water 


(E) The total or titration alkalinity (TA), which can be seen as the excess of proton 
acceptors over donors of a solution and is normally derived from an acidi- 
metric titration. The total alkalinity of a solution is relatively easy to measure, 
and which has the pleasant characteristic that it behaves conservatively as two 
waters are mixed or if temperature or pressure changes, similar to DIC. 
A formal definition and detailed treatment follow below. 


Let us revisit the simple CO2-H2O system from above with five species 
({H>CO3], [HCO; ], [CO3” ], [OH ] and [H*]). Titration of this system with a 
strong acid results in a decrease in pH and transformation of carbonate to bicar- 
bonate and bicarbonate to carbonic acid. Each bicarbonate ion can accept one 
proton, each carbonate ion consumes two protons before it is transformed into 
carbonic acid, and carbonic acid cannot accept protons. For this system, we can 
define HCO; as our reference (zero level of protons). Similarly, the hydroxide ions 
in solution will accept one proton, while any proton available just adds to the final 
proton concentration at the end of a titration. TA was formally defined by Dickson 
(1981) as the excess of proton acceptors ([HCO3 ] + 2 [C03] + [OH ]) over 
proton donors [H*] of the initial solution: 


TA = [HCO;] + 2[CO}"] + [OH] — [H+] (5.11) 
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The reader might have noted that this TA definition has much similarity to the 
charge balance and proton balance Eqs. (5.9), (5.10a, b) for this simple CO,—H,0 
system. If fact by re-arranging the charge balance (5.9) we obtain the charge balance 
alkalinity or excess negative charge (ENC): 


ENC = [HCO;] + 2[C0}] + [O0H7] — [H+] (5.12) 


and we see that this excess negative charge is equivalent to the titration alkalinity 
for this system. This excess negative charge introduced by Soetaert et al. (2007) is 
also known as the explicit conservative equation of total alkalinity (Zeebe and 
Wolf-Gladrow 2001; Wolf-Gladrow et al. 2007). 

Although seawater is a much more complex system, because it contains many 
potential proton acceptors and donors, such as borate, phosphate and silicate spe- 
cies, the same approach can be used. Dickson (1981) formally defined total alka- 
linity of filtered seawater by the excess of proton acceptors over proton donors with 
respect to the proton condition at pH 4.5 (Dickson 1981), which then leads to the 
expression: 


TA =[HCO;] + 2[CO;"] + [B(OH),] + [OH] + [HPO; ] + 2[PO;"] + [H3Si0;] 
+ 2[H)SiO7 | + [NH3] + [HS7] + 2[S?-] — [H+] — [HF] — [HSO;] 
—2[H2SO,] — [H3PO4] — [HNO] — [HNO3] 


(5.13) 


Following Wolf-Gladrow et al. (2007) and Soetaert et al. (2007), some additional 
proton acceptors (H>SiO,” , s7) and proton donors (HNO;, HNO3, H2SO,) have 
been added to the original definition of Dickson (1981). Most of these additional 
terms are close to zero in oxic seawater. Moreover, TA in marine waters is dom- 
inated by the carbonate species and simpler relations than (5.13) are often used. The 
carbonate alkalinity (CA) is defined as the alkalinity contribution of bicarbonate 
and carbonate species: 


CA = [HCO;] + 2[COz] (5.14) 


and usually accounts for >95% of the alkalinity, with borate contributing 
another ~3-4%. While the presence of these additional species hardly impacts 
measurement of TA by titration, it will affect the use of TA in calculations to obtain 
the sought-after distribution of bicarbonate and carbonate concentrations. 

The excess negative charge concept, as introduced earlier (5.12) can also be 
applied to seawater (Soetaert et al. 2007; Wolf-Gladrow et al. 2007): 
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ENC =[HCO;] + 2[CO}"] + [B(OH),] + [OH] + [H)PO,] + 2[HPO}] + 3[PO} ] + [H3SiO;] 
+ 2[H2SiOz-] + [NO;] + [No3] + [F7] + [HS7] + 2[S*”] + [HSO;] + 2[$07] — [H*]— [NH;]. 
(5.15) 


ENC and TA are related via: 


TA = ENC + X` NHy— $_ NO; - X_ NO, - X PO, -2X S0,- SOF 


(5.16) 


where the various -terms refer to total concentration of the respective acid-base 
pairs (e.g. JONH; = [NH,*] + [NH3]). This difference between TA and ENC is 
caused by components for which the species used as zero proton levels are charged. 
In other words, DIC and 5~B(OH); do not appear because they are uncharged at 
their zero proton level. This distinction between TA measured by titration (sensu 
Dickson) and ENC (charge balance alkalinity) is needed to understand the impact of 
biogeochemical processes on pH or pCOs, in particular those involving 
ion-exchange or nutrient uptake (Zeebe and Wolf-Gladrow 2001; Soetaert et al. 
2007; Wolf-Gladrow et al. 2007). For instance, nitrate or ammonium uptake by 
phytoplankton has a distinct effect on TA (decrease and increase, respectively; 
Brewer and Goldman 1976; Goldman and Brewer 1980) that cannot be inferred 
from the Dickson definition of TA (5.13), but is clear from Eqs. 5.15 and 5.16. The 
concept of potential alkalinity (Brewer et al. 1975) and its derivatives, such as TA*, 
are special cases of the ENC concept. The need to correct for nutrient cycling when 
using the alkalinity anomaly technique to quantify calcification is also made explicit 
when using Eq. 5.16. 


5.4 Buffering 


Natural waters, in particular seawater, contain multiple weak acids and their con- 
jugate bases. Any addition or removal of an acid or base results in re-adjustment of 
the equilibrium distribution of acids and conjugated bases, with the consequence 
that the disturbance is attenuated (the principle of le Chatelier). Buffering can 
involve only dissolved components (homogenous) or also solid phases (heteroge- 
neous). Heterogenous buffering in the ocean is called ocean carbonate compensa- 
tion (Box 5.2). Homogenous buffering of waters is most often quantified by the 
buffering value (B), which is formally defined as the concentration of acid or base 
to be added to influence pH (Van Slycke 1922; Urbanksy and Schock 2000): 


dC dC, 
Soe, (5.17) 
dpH dpH 
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where C, and C, are the concentration of added base and acid, respectively. The 
term buffer value as used here is identical to the terms buffer capacity, intensity and 
index in use by analytical chemists, geochemists and engineers. Buffer values are 
always positive because solutions resist changes according to the le Chatelier 
principle; consequently, a minus is needed when considering acid addition. The 
buffer value is a continuous non-linear function, as it is the inverse of the slope of 
the titration curve, i.e. a derivative (Stumm and Morgan 1996). The buffering value 
of seawater is expressed using the total alkalinity: 


rey (5.18) 
dpH 

For notational simplicity, we present buffer value as a derivative, while it is a partial 

derivative because all other concentrations are kept constant when determining this 

buffer value by titration (Hagens and Middelburg 2016). Van Slycke’s original 

definition of buffer value was in pH, but in seawater it is instructive to express 

buffering value in terms of proton concentration (): 


dTA 


B= -r (5.19) 


Buffer values B and ß are linked via B = —-In(10). H* « B. Table 5.2 gives the 
relevant equations for B (in pH) and ß (in proton concentrations). 

Total alkalinity includes multiple acid-base systems, and the overall buffering 
value can be obtained by summing the contributions of all the acid-base systems 
involved (Urbanksy and Schock 2000). The buffer value of seawater is mainly 
determined by the borate and carbonate contributions, and we therefore limited our 
calculations to these (and of course the contribution of water). The buffer value of 
seawater (S = 35, T = 15 °C, DIC = 2000 uM) is a function of pH and has two 
distinct maxima at pH 6 and 8.9, near the pK, and pK, of the carbonate system in 
seawater and a minimum at pH 7.4 (approximately midway between these maxima; 
Fig. 5.3). Within the pH range of 7.5-8.5 the buffer value increases with pH and 
can be approximated by the relation: 


K»DIC 
H+’ 


B = In(10) (5.20) 
showing that the buffer value of seawater in its normal range primarily depends on 
the DIC concentration, the pH and pK, value. Borate contributes about 15-22% to 
seawater buffering in this range. 
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Fig. 5.3 Buffer capacity of seawater. Top panel: total capacity, carbonate, borate and water 
contributions. Lower panel: total capacity, carbonate and borate contributions and approximate 
buffer capacity (see text) 


Another useful term related to the buffering is the sensitivity factor defined as the 
inverse of the buffer value expressed in terms of pH or protons (Soetaert et al. 2007; 
Egleston et al. 2010; Hagens and Middelburg 2016): 


dpH zj 
=B 21 
aTA a 
dH* 1 
— =p. 22 
ma? (5.22) 


These sensitivities and their use will be discussed below. 


5.6 Dissolved Inorganic Carbon Systematics 89 
5.5 Carbonate Mineral Equilibria 


In the above section, we only considered homogenous buffering, i.e. the 
re-equilibration among dissolved species. In natural systems, solid phases also 
contribute to buffering via dissolution and precipitation, and we call this hetero- 
geneous buffering. Carbonate minerals play a dominant role in heterogeneous 
buffering because of their omnipresence, high solubility and high reactivity. 

The precipitation and dissolution of calcium carbonate is represented by the 
reaction: 


Cat + CO} = CaCO; (5.23) 
for which we can write the equilibrium relation: 
Ky = [Ca ] [C03 | (5.24) 


where Ksp is the stoichiometric solubility product for calcium carbonate, which is a 
function of temperature, salinity, pressure and the calcium carbonate mineral. The 
solubility product is used to calculate the saturation state of seawater (Q) with 
respect to the calcium carbonate mineral: 


Q= bal ca (5.25) 


If Q = 1, then the solution is in equilibrium with that mineral, if Q > 1, then the 
solution is supersaturated with respect to that mineral and if Q < 1, then the 
solution is undersaturated with respect to that mineral. Undersaturation and 
supersaturation promote dissolution and precipitation, respectively. 

The two dominant calcium carbonate minerals in the ocean are calcite and 
aragonite, which differ in crystal structure, mineralogical properties (e.g. density) 
and solubility (Morse et al. 2007). The thermodynamic solubility products of calcite 
(pK = 8.48) and aragonite (pK = 8.30) differ by a factor 1.5, while their stoi- 
chiometric solubility products at salinity 35 and T = 25 °C are more than two 
orders of magnitude higher (pK = 6.36 and 6.18, respectively). The crystal struc- 
ture of carbonates has consequences for the incorporation of other ions, with larger 
cations (e.g., Sr°*, Ba’*) preferably incorporated in the orthorhombic aragonite 
structure and smaller cations (e.g., Mg?*, Mn**) fitting better in the hexagonal 
calcite structure. Dissolved magnesium concentrations in seawater are rather high 
(>5 times Ca?*) and Mg incorporation is therefore substantial (up to ~25% for 
biogenic carbonates), with large consequences for solubility. The solubility of 
Mg-calcite with about 2% Mg is lower than that of pure calcite, while calcites with 
about 12-15% Mg substitution are five times more soluble than aragonite 
(Arvidson and Morse 2014). 
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5.6 Dissolved Inorganic Carbon Systematics 


The conservative properties TA and DIC are often used as master variables to 
analyse the impact of environmental change on the marine CO, system, because the 
distribution of the species carbonic acid, bicarbonate and carbonate, and pCO, and 
pH depend on temperature and salinity. Figure 5.4 shows the temperature depen- 
dence of pH, pCO», bicarbonate and carbonate ion concentrations for seawater 
(TA = 2300 uM, DIC = 2000 uM). Increasing temperatures will cause an almost 
linear decrease in pH with a slope of —0.0155 pH/°C, while it will cause an 
exponential increase in pCO, values. Bicarbonate and carbonate ion concentrations 
will decrease and increase with increasing temperatures (slopes of about —0.55 
and +0.42 uM/°C, respectively). 

There are also distinct dependencies with salinity (Fig. 5.4), but these only include 
the impact of salinity on stoichiometric equilibrium constants, not the larger effect of 
salinity on total alkalinity via charge balance. Increasing salinity results in a decline in 
pH and carbonate ion concentration and an increase in bicarbonate and pCO>. 

Relations among TA, DIC, pCO, and pH are often presented graphically in the 
form of TA versus DIC plots with isolines for pH or pCO, (Fig. 5.5) for seawater 
(T = 15, S = 35). The addition or removal of TA or DIC or any of its constituents 
on such a plot is a vector property (Deffeyes 1965). For instance, CO, addition will 
increase DIC concentration, but not affect TA and thus result in a horizontal vector 
pointing to the right; pCO, of the water will increase and the pH will decline. 
Addition of protons will decrease TA, but not impact DIC; the downward heading 
vertical vector indicates that pH will decline and pCO, will increase. The addition 
of HCO; will impact both DIC and TA by one unit; the vector will have a slope of 
one and both pH and pCO, will slightly increase. The addition of CO,” will 
increase DIC by one unit and TA by two units, and the resulting vector indicates a 
lowering in pCO, and increase in pH. 


5.7 The Impact of Biogeochemical Processes 


Most biological and physical exchange processes occur on longer time scales than 
the chemical equilibrium reactions, and any addition or removal of a compound due 
to biology and physics thus causes re-equilibrations following the le Chatelier 
principle; i.e. buffering occurs. The impact of biogeochemical processes on the 
dissolved inorganic carbon system can be analysed (graphically or numerically) 
using their impact on DIC and TA. 

Gas exchange. Invasion of CO, from the atmosphere into seawater causes an 
increase in DIC, but does not change TA, with the consequence that pCO) increases 
and pH declines. The reverse process, an efflux of CO, from seawater to the 
atmosphere, results in pCO, decline and pH increase (Fig. 5.6). These CO2 
exchange processes may be initiated either by changing atmospheric mixing ratios 
of CO, or by cooling or warming of seawater. 
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Fig. 5.4 a Dependence of seawater pCO, and pH on temperature; b Dependence of carbonate 
and bicarbonate ions on temperature; c Dependence of seawater pCO, and pH on salinity; 
d Dependence of seawater carbonate and bicarbonate ions on salinity. Dashed lines are regression 


equations 
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Fig. 5.5 TA-DIC plot with pH contours (a) and pCO, (uatm) contours (b). Vectors are given for 
addition of protons, carbon dioxide, bicarbonate and carbonate 
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Fig. 5.6 Vector diagram on TA-DIC plot showing pH changes due CO, invasion and effluxes, 
carbonate mineral dissolution and precipitation, and primary production using nitrate as nitrogen 
source 


Carbonate minerals. Precipitation of calcium carbonate via the reaction: 
Ca** + 2 HCO; > CaCO; + CO, + H20 (5.26) 


removes one unit DIC from solution (-2 HCO, + CO2) and two units of TA (—2 
HCO; ), with the result that pCO, increases and pH declines. Dissolution of cal- 
cium carbonates results in the release of one unit DIC and two units TA and thus a 
decline in pCO, and pH increase (Fig. 5.6). 

The above examples directly involved one of the species of the dissolved 
inorganic carbon system, but many important biological processes do not directly 
involve any of these (e.g., oxidation-reduction reactions, such as nitrification or 
sulphide oxidation) or also involve other compounds besides the inorganic carbon 
species (e.g., primary production and organic matter degradation), and these do 
impact pH, pCO, and other parameters. To include these additional processes, we 
have to return to the relation between Excess Negative Charge (ENC) and Titration 
Alkalinity (TA). 
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Primary production involves the fixation of inorganic carbon and the assimila- 
tion of nitrogen and phosphorus to produce Redfield organic matter (see Chap. 6). 
Nitrogen assimilated in the form of either nitrate or ammonium has consequences 
for alkalinity because of the principle of nutrient-proton compensation mechanism 
(Brewer and Goldman 1976; Soetaert et al. 2007; Wolf-Gladrow et al. 2007). The 
rationale is as follows. Uptake of a nutrient in ionic form requires either 
uptake/release of a proton or hydroxide ion if internal compensation (e.g. NH4* 
with PO,° ) is not occurring, otherwise the organism would be charged. Ammo- 
nium uptake results in a decrease in alkalinity because either a proton is released or 
an additional hydroxide ion is assimilated for charge compensation. By analogy, 
nitrate uptake increases alkalinity because either a proton is consumed or a 
hydroxide ion is released to maintain electroneutrality. 

Accordingly, primary production based on ammonium can be represented as: 


CO, + n NH; + p H3;PO4+ H20 => (CH20) (NH3), (H3PO14), + © (5.27) 


where n and p are the stoichiometric ratios (16/106 and 1/106 for Redfield ratios, 
respectively). This reaction lowers TA by n-p, i.e. 15/106 equivalents and DIC by 
one unit. 

Similarly, primary production based on nitrate, i.e., 


CO2 + n HNO; + p H,;PO4+ (1 +n) H20 
= (CH20)(NH3)„(H3PO14), + (1 + 2n) O2 (5.28) 


increases TA with n + p, i.e. 17/106 equivalents and DIC by one unit (based on 
Redfield). 

Note that these reactions have been presented in the form of CO2, NH3, HNO3 
and H3PO, for convenience, but could have been written also in terms of bicar- 
bonate, ammonium, phosphate and nitrate ions. Because of the electroneutrality 
condition and rapid equilibration, it does not matter which species is actually taken 
up. For instance, when CO, is taken up during photosynthesis, buffering will 
replenish CO, from the large bicarbonate pool, thereby consuming protons, i.e. 
increasing the pH. Alternatively, when bicarbonate is assimilated, electroneutrality 
maintenance would imply uptake of a proton or release of a hydroxide ion, the 
result being an increase in pH. 

As another example, consider the anaerobic oxidation of methane coupled to 
sulphate reduction. We can write this equation either as 


CH; + H,SO4 > H2S + CO, + 2 H,O (5.29a) 
or 


CH4 + SO% => HS7 + HCO; + H20 (5.29b) 
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Table 5.3 Biogeochemical processes and changes in excess negative charge and total alkalinity 
(Soetaert et al. 2007) 


Process Reaction AENC ATA 
Aerobic mineralization (CH20) (NH), (H3POs), +02 < CO2 + n NH; + p H;PO4 + H20 0 n-p 
Denitrification (CH20)(NH3),(H3PO,), + 0.8 HNO; = CO) +n NH; + p H,PO, +0.4N3 + 0 0.84 n— 
1.41H,0 p 
Mn-oxide reduction (CH20) (NH3), (H3PO4), +2 MnO, + 4H* <> CO) + n NH; + p HPO, + 4 n-p +4 


2 Mn’+ + 3H,0 


Fe-oxide reduction (CH,0)(NH3),, (H3POx), +2 Fe20; + 8H* & CO, + n NH; + p H,PO4 + 8 n-p +8 
4 Fe’* + 5M0 


Sulfate reduction (CH20) (NH), (H3PO4), +0.5 H2SO0, < CO, + n NH; + p H3PO4 + 0 n-p+1 
0.5 H2S + H20 

Fermentation (CH20)(NHG),,(H3PO4), < 0.5 CO2 + n NH; + p HPO, + 0.5 CH, + H20 0 n-p 

Anaerobic oxidation of CH, + H2S0, = CO: + H2S + 2 HO 0 2 

methane 

Calcite precipitation Ca? + CO = CaCO; 2 2 

Primary production CO + n HNO; + p H;PO, + (1 + n) H20 + (CH20)(NH;3),(H3PO4), + 0 p+n 

(nitrate) (1 + 2n) © 

Primary production CO + n NH; + p HPO, + H20 + (CH20)(NH;)„(H3PO4), +02 0 p-n 


(ammonium) 


n = N/C ratio of organic matter and p = P/C ratio of organic matter 


Irrespective of the way formulated TA increases by two units and DIC by one 
unit. Soetaert et al. (2007) have worked out in detail the impact of biogeochemical 
processes on changes in ENC and TA (Table 5.3) and the related pH change. These 
pH changes are dependent on the actual pH of the system. 

To illustrate the latter, we return to the precipitation of calcium carbonate. The 
precipitation of calcium carbonate can be presented as: 


Ca?’ t + CO3 => CaCOs, (5.30a) 

Ca?t + 2 HCO; = CaCO3 + CO: + H20, (5.30b) 
Ca? t + HCO; => CaCO; + Ht (5.30c) 
Ca** + HCO; => CaCO; + 2H. (5.30d) 


While all these reactions are correct in terms of mass balance and stoichiometry: i.e. 
two units TA and one unit DIC are consumed for precipitation of one mole CaCOs, 
the second one indicates production of CO, and the last two imply production of 
protons. This non-uniqueness is the result of re-equilibration reactions, and these 
have to be taken into account when presenting calcium carbonate precipitation or 


96 5 Biogeochemical Processes and Inorganic Carbon Dynamics 


dissolution in a single equation. Evidently, at pH values >pK3, the carbonate ion is 
the most important species (Fig. 5.1) and reaction (5.30a) will be the dominant 
reaction. This is the equation most often used in laboratory studies performed at 
high pH and no protons or carbon dioxide are net generated. Most field studies 
employ reaction scheme (5.30b) and have shown that calcium carbonate precipi- 
tation represents a source of carbon dioxide, but not one mole carbon dioxide for 
one mole of calcium carbonate because of buffering. Frankignoulle et al. (1994) 
presented the factor which expresses the amount of carbon dioxide generated per 
unit calcium carbonate precipitated: 


Ca** + (1+) HCO; > CaCO; + WCO,+ (1 — y) HY. (5.31) 


The factor y basically merges reactions (5.30b, c) above and can be calculated 
analytically and has a value of 0.6—0.7 for most marine waters. The result that pCO2 
increases and pH declines is consistent with the vector analysis presented above 
(Fig. 5.5). 

An alternative generic treatment has been presented by Hofman et al. (2010). 


Ca?* + agH»CO3 + a; HCO; +a,CO} = CaCO; + (2a +a) H} (5.32) 


where o = H,CO,/DIC, a, = HCO3 /DIC, a = CO, /DIC and op + o1 + 
%2 = 1. The number of protons released by calcite precipitation is given by 20 +41, 
i.e. two times the carbonic acid + bicarbonate contribution to the DIC pool and 
varies from 2 at low pH to zero at high pH values (Fig. 5.7, red dashed line). 
However, the protons generated during calcite precipitation will be buffered, and 
the actual increase in proton concentration is given by eee i (Hofman et al. 


2010; red solid line) where B is the buffer value in terms of protons introduced 
earlier (Eq. 5.22). The released protons will decrease the pH of the solution (blue 
dashed line). A similar analysis can be made for calcite dissolution: two protons are 
consumed during calcite dissolution at low pH, the net proton change is less 
because of buffering and the pH increase shows a non-linear response with two 
maxima (blue solid line). 

Using the information in Table 5.3 and the approach of Soetaert et al. (2007), it 
is possible to calculate the impact of any process on pH and how it varies with pH. 
The pH dependence is the product of the change in negative charge (i.e. number of 
protons involved) and the sensitivity factor (Eq. 5.21). Figure 5.8 shows the pH 
dependence of the sensitivity factor with distinct maxima at pH 4.3, 7.3 and 9.9. 
Aerobic mineralization releases carbon dioxide, ammonium and phosphate. Carbon 
dioxide production decreases the pH for pH > ~5 (Fig. 5.9), in particular at pH 9.9 
and 7.3 because of poor buffering, but pH increases for pH < ~5 because of 
ammonium release. The release of carbon dioxide does not impact the negative 
charge (or proton balance) at low pH because all inorganic carbon is already present 
as carbon dioxide. Denitrification lowers pH for pH > ~7, while it increases pH for 
pH < ~7 because of nitrate consumption and ammonium release, both impacting 
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Fig. 5.7 Calcium carbonate precipitation/dissolution, proton and pH changes as function of pH. 
The number of proton generated per mole calcite formed/dissolved before buffering (red dashed 
line) and after buffering (red solid line). The decline in pH for precipitation of one umole of calcite 
(dashed blue line) and increase in pH for dissolution of one umole calcite (solid blue line). There is 
a break at the pH where calcite starts to dissolve or precipitate 


the negative excess charge and TA (Eq. 5.16), in particular at pH 4.3 because of the 
high sensitivity factor (Fig. 5.8). Manganese and iron oxide reduction always 
increase pH with the relative pH increase primarily depending on the sensitivity 
factor (Fig. 5.9). 


Box 5.1: Ocean acidification 

Atmospheric carbon dioxide concentrations steadily increased during the last 
century because of the use of carbon-based energy resources, changes in land 
use, and lime production. Part of this additional carbon dioxide remains in the 
atmosphere, i.e. an airborne fraction of about 45%, and is the main driver of 
global warming. The other part ends up in the ocean or terrestrial biosphere. 
The ocean uptake of anthropogenic carbon accounts for 25-30%, but this 
service to humankind comes at a price: acid-base equilibria in the ocean have 
shifted. DIC and bicarbonate concentrations increased (carbonation), while 
carbonate ions and pH declined (ocean acidification). Ocean acidification, or 
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the “other CO, problem”, has consequences for the chemistry, biology and 
geology of the ocean. 

Figure 5.10 shows the dependence of pH and DIC on atmospheric pCO, 
for a constant alkalinity ocean (2300 uM, T = 15, S = 35). Dissolved inor- 
ganic carbon concentrations increase almost linear with a slope of 0.513 uM/ 
patm pCO, (or 13.745 uM/uM dissolved CO3), while pH declines almost 
linearly with a slope of —0.0011 pH/uatm pCO, (or —0.027 pH/uM dissolved 
CO2). These thermodynamic predictions for declining pH and increasing DIC 
are fully consistent with observations in the ocean. These dependences of 
CO, and pH have received much attention and have been treated more for- 
mally using differential calculus (Sarmiento and Gruber 2006; Hagens and 
Middelburg 2016). The changes in pH can be split into multiple parts: 


_ (pH _ (pH _ (OpH _ (pH 
dpH = (Fe) ar (Fees. (Fe) ara Gan dDIC 4 


(5.33) 
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Fig. 5.9 a pH change due to iron-oxide and manganese oxide reduction as function of pH; b pH 
change due to aerobic mineralization and denitrification as a function of pH. The sensitivity factor 
(Fig. 5.8) is shown as reference (grey) 


where (4) are partial derivatives, implying that these are, in principle, only 


valid if the other variables are kept constant, and they can be considered as a 
measure of the sensitivity of pH to a change in the respective environmental 


variable. For instance, the first partial derivative, (2) has a value of 
~-0.014 per degree (Hagens and Middelburg 2016), close to slope in 
Fig. 5.10 (-0.015 per degree). The third one, (#4) , is the sensitivity factor 
(Eq. 5.21) shown in Fig. 5.8. 
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Fig. 5.10 pH and DIC dependence on pCO, (1atm) (a) and dissolved CO, (uM) (b). Regression 
lines are shown as dashed lines 


Ocean acidification will lead to a decrease in buffer value (Fig. 5.3), and 
this is one of the reasons why the ocean is projected to take up less 
anthropogenic carbon in the future (besides changes in ocean physics and 
biology). Consequently, the ocean carbon dioxide system will be more sen- 
sitive to changes, and diurnal and seasonal changes in pH are projected to 
increase. Moreover, shifting chemical equilibria will have consequences for 
sound attenuation in the ocean. 
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Many organisms, autotrophs as well as heterotrophs, will be impacted by 
either the decline in pH and carbonate ions or the increase in dissolved 
inorganic carbon. The latter, carbonation, might stimulate primary producers, 
which are carbon limited at the moment. Most calcifying organisms are likely 
going to suffer from the decrease in carbonate ion availability, but calcifi- 
cation is under biological control and sometimes completed inside the 
organism, such that these responses are rather complex (Kroeker et al. 2013). 


Box 5.2: Carbonate compensation dynamics 

Buffering in the ocean not only occurs homogenously by re-arrangement of 
acids and conjugated bases, but also heterogeneously by re-adjustment of the 
balance between precipitation and dissolution of carbonate minerals. While 
homogenous buffering is fast (instantaneous equilibria), heterogeneous 
buffering has a larger capacity because of the large stock of carbonate min- 
erals stored in marine sediments. This heterogeneous buffering is called 
carbonate compensation and biology plays a major role. 

Carbonate minerals are predominantly formed biologically in the modern 
ocean, but there are a few exceptions, such as the formation of ooids in 
tropical systems and authigenic calcite and dolomite formation in sediments. 
However, even in these cases, biology plays indirectly a major role by 
governing the chemical composition of the fluids they formed in (e.g. car- 
bonate formation induced by alkalinity production resulting from anaerobic 
oxidation of methane). Biological carbonate formation takes place in the 
water column by autotrophs (e.g., coccoliths) and heterotrophs (pteropods 
and foraminifera) and in the benthos by various organisms, autotrophs (e.g. 
coralline algae) and heterotrophs (e.g., corals, crustaceans, and molluscs). 
Biogenic carbonates can be aragonite, calcite or high-Mg calcites and com- 
binations thereof. 

Following death of calcifiers in benthic systems, the biogenic carbonate 
can either dissolve or accumulate at the seafloor and be buried. Carbonate 
produced in the surface layer of the open ocean can dissolve in the water 
column while particles settle, dissolve at the seafloor, or accumulate in sed- 
iments. Accumulation of biogenic carbonate is a prominent feature of ocean 
sediments; some sediments consists almost entirely of biogenic carbonate 
debris. Dissolution of carbonate is primarily driven by undersaturation, and 
organisms contribute to dissolution in a number of ways. One, many biogenic 
minerals have organic layers and microbial degradation of these layers 
exposes new surfaces to undersaturated solutions. Two, boring organisms 
(sponges, fungi) weaken the structure and texture of biogenic carbonate and 
consequently accelerate dissolution. Three, the metabolic activity of organ- 
isms has consequences for the saturation state of solution, e.g. carbon dioxide 
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Fig. 5.11 Carbonate compensation concept. Biogenic CaCO; tests are produced in the photic zone 
of the oceans (green circles). Upon death, those tests escaping dissolution near the surface, settle 
along with clays materials. Above the saturation horizon, waters are supersaturated and CaCO; tests 
are largely preserved. Below the saturation, waters are undersaturated because of increasing 
solubility with depth and the release of CO, from organic matter decay and CaCO; will dissolve. 
Dissolution occurs primarily at the sediment surface as the sinking velocity of debris is rapid (broad 
white arrows). At the carbonate compensation depth, the rate of dissolution exactly matches rate of 
supply of CaCO; from above. At steady state this carbonate compensation depth is similar to the 
snowline; the first depth where carbonate poor sediments occur. The lysocline is the depth interval 
between the saturation and carbonate compensation depth (from Boudreau et al. 2018) 
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release during respiration and generation of strong acids during re-oxidation 
reactions (nitrification, sulphur oxidation) lowers the saturation state. The 
latter processes are important in coastal sediments that receive high organic 
carbon inputs and that are bioturbated or vegetated. 

There are two carbonate compensation mechanisms operating in the 
ocean: chemical and biological compensation (Boudreau et al. 2018). 
Chemical compensation focuses on the dissolution or preservation of car- 
bonate minerals at the ocean floor and implicitly assumes that net carbonate 
production remains constant. Carbonate particles settling in the ocean interior 
start to dissolve when water becomes thermodynamically undersaturated at 
depth Zsa. This carbonate saturation depth can be estimated from: 


[Ca?+ ] cor 


Zsat = Zrefln ( K (5.34) 
sp 
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where K,, is the temperature, salinity and pressure depending stoichiometric 
solubility product, [Ca] and [C07] are the concentration of dissolved 
calcium and carbonate ions and Z,ef is a scaling parameter (Boudreau et al. 
2010). This undersaturation is partly due to increasing pressure and declining 
temperature (thermodynamics) and partly the consequence of a decrease in 
carbonate ion because of respiration of particulate organic matter (the bio- 
logical pump). Below z,,, dissolution rates increase systematically with depth, 
and at a certain depth the dissolution rate balances the settling flux of car- 
bonate (Fig. 5.11). This carbonate compensation depth (Zccp) is governed by 
the following equation (Boudreau et al. 2010): 


(5.35) 


Fear[Ca°*] | EAN 
KopAP mi Ksp 


Zccp © castn( 


where Fear is the export flux of carbonate, A is the surface area of the seafloor 
and Bm is the mass transfer of solutes across the diffusive boundary layer at 
the seafloor. The carbonate compensation depth (Zccp) is always larger than 
the saturation horizon (Zsa) because of the first term with only positive 
parameters. The carbonate compensation depth is equal to the snowline, the 
depth at which carbonate disappear from sediments, under steady-state con- 
ditions. During periods of ocean acidification, bottom-waters will eventually 
obtain lower carbonate ions concentrations, and the saturation and carbonate 
compensation depths rise and dissolution of carbonate at the seafloor 
increases until a new balance between dissolution and export flux of car- 
bonate has been reached. Conversely, during periods of alkalinisation, bottom 
water will eventually get higher carbonate ion concentrations (by advection), 
carbonate mineral dissolution decreases, resulting in a deepening of saturation 
and compensation depth till a new balance has been reached. This carbonate 
compensation mechanisms operates on a time scale of 100—10,000 years. 

This chemical compensation mechanism is based on the assumption that 
ocean acidification has no impact on calcification and the export of calcium 
carbonate. However, many experiments have shown that calcification rates, 
and thus export of calcium carbonate from the surface ocean declines with 
saturation state of surface water. A decline in carbonate export (Fear) would 
lead to shallowing of the carbonate compensation depth on the very short 
time scale (<1 yr), but would cause additional deepening on the longer term 
(>10* yr) because calcification is an alkalinity sink. With less removal of 
carbonate ions in the surface waters, deep water will eventually become richer 
in carbonate and more carbonate minerals will survive dissolution at the 
seafloor (Boudreau et al. 2018). 
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Organic Matter is more than CHO 


In the previous chapters, we have presented organic carbon in isolation from the 
other elements and focused on the biological processes involved in the transfor- 
mation from inorganic to organic carbon and vice versa. This approach neglects the 
fact that the element carbon is part of molecules. The thermodynamic stability and 
reactivity of, and interactions among, these molecules eventually determine the 
function of organic carbon compounds and govern the rates of organic matter 
production, transformation and consumption. Organic molecules contain not only 
C, H and O, but also N, P, S and multiple other elements, and often in certain ratios, 
e.g. one P atom in each DNA or RNA nucleotide. Organic molecules have func- 
tional groups and stereochemistry that determine their interactions with the 
environment. 

Following a concise discussion of Redfield organic matter and non-Redfield 
organic matter, we focus on organic matter as food for organisms, the composi- 
tional consequences of preferential consumption and the consequences for the 
composition of organic matter preserved and buried in marine sediments. 


6.1 Redfield Organic Matter 


Chemical analysis of marine particulate organic matter, organisms and seawater by 
early oceanographers inspired Redfield (1934) to propose that marine plankton have 
relatively constrained atomic ratios of ~140 C to ~20 N and ~1 P (Fig. 6.1). 
This empirical Redfield ratio was originally based on the similarity of these ratios in 
seawater and in marine plankton, and these ratios have since been shown to apply to 
assimilation ratios during phytoplankton growth and regeneration ratios during 
marine organic matter degradation. Redfield further noted the implications for 
oxygen consumption and later on (Redfield 1958) posed the intriguing question 
whether the biological processes control the proportions of these elements in the 
water. 
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Fig. 6.1 Redfield’s original 30 
correlation between 

phosphate and nitrate in 
seawater (Redfield 1934). The 25 
line has a slope of 20 rather 
than the modern value of 16 
for ideal Redfield 20 


The Redfield ratio as formalized by Redfield et al. (1963), at 106C:16 N:1P, has 
been validated by thousands of observations, and the overall ratios appear to be 
robust. Redfield ratios have shown useful in a wide range of topics, from nutrient 
limitation to biogeochemical modelling in the present and past ocean. They are used 
in reconstruction of anthropogenic carbon inventories and are in the equations 
underlying Earth System modelling (Sarmiento and Gruber 2006). The original 
Redfield ratio has been extended to include other elements and heterotrophic 
organisms (ecological stoichiometry) and has been related to cellular level and 
global scale processes (Falkowski 2000; Sterner and Elser 2002; Box 6.1). At the 
cellular level, nitrogen is primarily used for the construction of proteins and 
phosphorus for the synthesis of ribosomal ribonucleic acids (rRNA). Loladze and 
Elser (2011) have shown that protein and ribosome synthesis results in a protein: 
tRNA ratio of ~3, which corresponds to an atomic N:P ratio of 16, as observed in 
the ocean for plankton. Moreover, they have shown that N and P limitation during 
cell growth results in N:P ratios below or above the Redfield ratio, respectively. 

Redfield et al. (1963) presented the (canonical) overall reaction for marine 
organic matter production: 


106 CO; + 16 HNO; + H3POq4 + 122 H20 = (CH20)j06(NH3) ;6H3PO4 + 138 O2 
(6.1) 


and consumption: 


(CH20)9¢(NH3) ,sH3PO4 + 138 O2 > 106 CO, + 16 HNO; + H3PO4 + 122 H20 
(6.2) 
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These RKR equations indicate that 138 units of oxygen are required for oxidation 
of one mole of Redfield organic matter: 106 oxygen molecules to convert CHO to 
CO, and 32 oxygen molecules to convert the NH; all the way to HNO3. The O2:C 
coefficient for aerobic respiration (including nitrification) is thus 1.3. Primary 
production based on nitrate (new production) or ammonium (recycled productions) 
has consequences for the quantity of oxygen released. 

The problem with the traditional Redfield ratios as formulated in Eqs. 6.1 and 6.2 
is that it simplifies marine organic matter into an elementary carbohydrate (CH,O), 
whereas newly produced marine organic matter is better represented as a mixture of 
proteins (~ 65%), lipids (~ 19%), and carbohydrates (~ 16%) and other compounds 
(pigments and nucleic acids; Hedges et al. 2002; Bianchi and Canuel 2011). While 
simple carbohydrates are the first molecules formed during carbon fixation, the 
biosynthesis of proteins, lipids and polysaccharides results in the loss of hydrogen 
and oxygen. Hydrogen and oxygen loss results from dehydration, wherein CHO is 
transformed into CH, chains in proteins and lipids. Moreover, there is also hydrogen 
loss during the transformation of the NH3-group in ideal Redfield organic matter to 
the NH-amide group in proteins. Estimates of marine organic matter based on 
mixing of proteins, lipid, carbohydrates and nucleic acids and using modern ana- 
lytical tools (Anderson 1995; Hedges et al. 2002), arrive at average compositions of 
C196H175042N16P (Anderson 1995) and C1096H177037N16PSo.4 (Hedges et al. 2002) 
rather than C1o6H2600106N16P of the ideal Redfield ratio. This lower oxygen content 
has consequences for the O,:C respiration coefficient. Equation 6.2 can be gener- 
alized to: 


C,HpO,N5P; + yO. = «CO, + OHNO; + £ HPO; + 2 H2O (6.3) 


where y =a + 0.25 B — 0.5% + 1.258 + 1.25 e and à = —0.5ò — 1.5 € + 0.5 B. 
Substitution of these revised marine organic matter compositions in Eq. 6.3 results 
in O2:P ratios of 150 to 154 and an O2:C respiration coefficient (y/a) of ~ 1.43. The 
latter result is consistent with the O:C ratio of ~ 1.46 obtained from inverse 
analysis of water-column nutrient and oxygen changes along isopycnals (Anderson 
and Sarmiento 1994). 


6.2 Non-redfield Organic Matter 


The composition of marine particles is rather uniform because phytoplankton and 
bacteria (and their remains) dominate plankton, and these organisms are primarily 
made up of proteins (>50%) supplemented by lipids, carbohydrates, nucleic acids 
and pigments. Individual species may deviate from ideal or revised Redfield ratios, 
but over larger spatial and temporal scales these differences average out. The 
composition of seagrasses and, in particular, emergent saltmarsh vegetation and 
mangrove trees deviates from unicellular primary producers: higher C:N and C:P 
ratios, less proteins and nucleic acids and more structural carbohydrates and lignins. 
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Moreover, there is much more variability in the composition of benthic primary 
producers, particularly in macrophytes. The higher variability is due to tissue dif- 
ferentiation: compounds involved in photosynthesis and carbon fixation are pri- 
marily in the leaves and not in the roots. Carbohydrates for storage (e.g. starch) can 
be transferred from leaves to rhizomes for long-term storage. These marine 
macrophytes re-allocate compounds and essential nutrients over the growing sea- 
son: e.g. the resorption of nitrogen from senescent leaves and transfer to active 
growing leaves. Structural carbohydrates and lignins are needed for strength to deal 
with currents, waves and wind. Lignins are cross-linked phenolic polymers that 
occupy the space between structural carbohydrates, such as cellulose in the cell wall 
of macrophytes and trees. Polysaccharides are permeable to water, while the more 
hydrophobic lignins are not. Crosslinking between lignins and celluloses thus not 
only strengthens the cell wall but also provides capability to transfer water effi- 
ciently from roots to the leaves, where evapotranspiration occurs. The 
lignin-carbohydrate provided strength and stiffness are needed for saltmarsh plants 
and trees to stand and for submerged vegetation to resist currents and waves. These 
polymers are relatively rich in carbon and poor in nitrogen and phosphorus; this is 
the reason why seagrass, salt-marsh plants and in particular mangroves have high 
C:N and C:P ratios. These polymeric substances are difficult to degrade, with the 
consequences that the palatability of seagrasses, marsh vegetation and mangroves is 
relatively low (see Fig. 3.10). 


6.3 Organic Matter is Food 


Organic matter, detritus as well as living biomass, is consumed because it is the 
primary source of energy and nutrition for heterotrophic consumers, i.e. it is food 
for animals and a substrate for heterotrophic microbes. Osmotrophic organisms 
either use dissolved organic matter directly or after extracellular (cell-attached or 
free) hydrolysis of polymeric substances into smaller units that can pass the cell 
membrane of these microbes. Animals ingest particulate organic matter, part of 
which is digested, and part is egested. The energy balance for animals is (Welch 
1968): 


I=G+R+E (6.4) 


where I is the ingestion of food, G is growth and reproduction of the animal, R is 
respiration and E is egestion and excretion (e.g. feces). Food ingested by animals is 
digested physically (e.g. chewing, grinding), chemically (by enzymes) and bio- 
logically (by microbes in the digestion organs) to make it amenable to enter tissues, 
where it can be used for growth and respiration (energy). This simple equation can 
be used to define assimilation (A = G + R), assimilation efficiency (A/I), gross 
growth efficiency (G/I) and net growth efficiency (G/A). 
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Moreover, during periods of sufficient food, consumers can internally store part 
of the food assimilated for reserves in the form of carbohydrates and lipids. 

Assimilation efficiencies vary widely depending on the food quality and the 
consumer: carnivores generally have higher assimilation efficiencies than herbi- 
vores and detritivores. However, this is partly compensated by higher net growth 
efficiencies for the herbivores, with the result that gross growth efficiencies vary 
over a rather narrow range (15-35%, Welch 1968). Assimilation efficiencies vary 
between elements, between biochemical classes of compounds (carbohydrates, 
amino acids, lipids) and among molecules within a class of compounds. For 
instance, carbohydrates and amino acids enriched in intracellular materials have 
higher assimilation efficiencies than those in cell-wall materials because cell-walls 
are less digested (Cowie and Hedges 1996). The composition of the diet and the 
material assimilated often differs from the composition of the consumers; conse- 
quently, some compounds (e.g. carbohydrates) are more used for respiration, while 
other compounds (e.g. amino acids) may be preferentially directed towards growth 
and synthesis of new tissues. 

Detailed feeding studies of lipids, carbohydrates and amino acids at the com- 
pound levels with marine animals have revealed that some compounds are pref- 
erentially consumed and used for respiration, while others are used to synthesize 
compounds de novo or by transformation of assimilated compounds. These patterns 
have been shown not only to be taxon specific, but also to have a microbial 
processing signature (Woulds et al. 2014). The latter is consistent with the 
microbiome concept, i.e. microbes within animals’ digestion system are key to the 
functioning of the animal. Moreover, there are a number of compounds that are 
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Fig. 6.2 Animal processing of organic matter. Food taken in is either digested or not and then 
egested. Digested organic matter is used for respiration or for growth that occurs primarily via 
synthesis of new compounds, but also via transformation or retention of compounds taken in. 
Moreover, some compounds can be produced de novo or transformed by microbes in the digestive 
system 
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considered essential because consumers cannot synthesize them and have to 
assimilate them from their diet or from the microbes in their digestive system. These 
include some poly-unsaturated fatty acids and amino acids, such as threonine, 
valine, leucine, isoleucine, phenylalanine and lysine. The combined effect of 
compound specific processing during digestion, the use of compounds for growth, 
storage or respiration, and the need to acquire essential compounds and elements 
results in changes in the organic matter remaining, excreted and produced by these 
consumers (Fig. 6.2). 


6.4 Compositional Changes During Organic Matter 
Degradation 


The above differences in assimilation efficiencies among various compounds and 
biochemical classes, and the preferential consumption of more labile components 
(Box 3.1; Fig. 3.9) have consequences for the composition of the organic matter 
remaining. Hedges, Wakeham, Lee and colleagues systematically studied the 
changes in the biochemical composition of organic matter during degradation from 
fresh phytoplankton, via sediment traps to surface and subsurface sediments 
(Wakeham et al. 1997). Amino acids (proteins), lipids and carbohydrates dominate 
the composition of phytoplankton and of detrital organic matter in shallow traps 
(Fig. 6.3). The composition changes rapidly upon degradation and organic matter in 
sediment traps below the surface mixed layer: amino acids decline from >60% to 
~ 25% and lipids from ~ 10 to ~2%, while the relative carbohydrate contribution 
remains similar because some carbohydrates are structural components. Moreover, 
the proportion of organic matter that cannot be characterized molecularly increases 
with progressive degradation from a few % in fresh phytodetritus to >50% in 
sediment trap organic matter and >70% in deep-sea sediment organic matter 
(Fig. 6.3). 

The majority of organic matter in deep-sea sediment traps and sediments cannot 
be characterized using solution or gas-based chromatographic methods because of 
low organic solvent extraction and hydrolysis yields. Solid-phase NMR techniques 
have been applied as an alternative, and these studies not only confirmed the 
proportions of amino acid and carbohydrate based on solution-based techniques but 
also revealed that non-hydrolysable carbon-rich material (i.e. black carbon) and 
non-protein alkyl group made up most of the molecularly uncharacterizable organic 
matter (Fig. 6.4). Accordingly, the organic matter buried in sediment not only 
represents a small fraction of that produced (few %, Chap. 4), but it also differs 
significantly and systematically from phytodetritus (Fig. 6.3). This complicates the 
use of bulk organic matter properties as a proxy for the origin of the organic matter 
and for reconstruction of paleoenvironments (Middelburg 2018). 

The susceptibility towards degradation varies systematically among biochemical 
classes: pigments > lipids = carbohydrate > amino acids > lignin > black carbon. 
Pigments, DNA, RNA and other cellular constituents are the most easily degradable 
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Fig. 6.3 Changes in the biochemical composition of organic matter during degradation from 
plankton, via sediment traps at different depths to surface and subsurface sediments in the Pacific 
Ocean (Wakeham et al. 1997). Uncharacterized fraction was calculated by difference: the 
uncharacterized fraction increases with progressive degradation 
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Fig. 6.4 Biochemical composition of marine sediments based on solid-phase NMR data (Gelinas 
et al. 2001). Four major compounds group were identified: amino acids, carbohydrates, black 
carbon and non-protein alkyl (AlkyINP) 
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compounds, followed by lipids, carbohydrate and amino acids (Fig. 6.3). While 
lignins are degraded in aerobic soils by fungi, these compounds are rather stable in 
the marine environments with the consequence that they have high preservation 
potential. Besides these differences among biochemical classes, there are also large 
differences among compounds within biochemical classes because of the inherent 
structural differences, as well as the macromolecular context. Structural carbohy- 
drates are more stable than those involved in storage. Amino acids incorporated in 
large proteins are less available for consumers than those dissolved as free amino 
acids. Lipid reactivity towards degradation varies widely, as it depends on multiple 
factors, including the structure, stereochemistry, head-group, degree of saturation, 
and ether or ester linkages. While the relative reactivity of various compound 
classes and compounds within a biochemical class is systematic, absolute rates of 
degradation are largely context depending, i.e. identical compounds may have order 
of magnitude differences in kinetic parameters depending on the environment. 
These systematic changes in organic matter composition due to degradation can 
be used to quantify the progress of degradation, or in other words, the degradation 
history. Often used degradation state parameters are intact to total pigment ratios, 
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Fig. 6.5 Systematic changes in organic matter composition as degradation progresses. The 
degradation index of Dauwe and Middelburg (1998) is based on subtle changes in the amino acid 
composition. Fresh phytoplankton material has positive values, while extensively degraded 
sedimentary organic matter in deep-sea sediments has values <—1.5. Degradation causes relative 
accumulation of a non-protein amino acids Bala and Gaba and b aminosugars and a decrease in 
c the organic carbon content per sediment surface area (OC/SA) and d the fraction present as 
non-protein alkyl-C. Data are from Dauwe et al. (1999, blue), Gelinas et al. (2001, red) and 
Vandewiele et al. (2009, black) 
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pigment to carbon ratios, the contribution of amino acid carbon or nitrogen to the 
total carbon or nitrogen pool, the fraction of organic matter present as alkyl-C and 
the relative contribution of carbohydrates and amino acids to the total carbon pool. 
The most generic degradation state proxy is the amino acid based degradation index 
of Dauwe and Middelburg (1998). This degradation index is based on subtle 
changes in the amino acid composition of particulate organic matter. The resulting 
degradation index varies from +1.5 for fresh phytoplankton material to —2.5 for 
intensively degraded deep-sea sediment organic matter. The rationale is that the 
amino acid composition of phytoplankton is rather uniform and that changes 
measured in particulate organic matter can be attributed to mineralization processes. 
The degradation index correlates with multiple other organic matter degradation 
proxies, such as accumulation of microbial detritus (aminosugars, bacterial degra- 
dation products, D-amino acids derived from bacterial cell walls) and the decrease 
in organic carbon per unit surface area and contribution of alkyl-C (Fig. 6.5). The 
degradation index provides a continuous parameter to infer the organic matter 
degradation history, and one would expect a relation with the first-order rate con- 
stant for organic matter degradation shown in Fig. 3.9. Progressive degradation of 
organic matter should, according to the reactive continuum concept, lead to lower 
reactivity rate constant and to a lower degradation index. Figure 6.6 shows that the 
logarithm of the first-order rate constant and the degradation index of bulk organic 
matter are indeed correlated and that first-order rate constant can, in principle, be 
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Fig. 6.6 Relation between the degradation index and the first-order reaction constant (data from 
Dauwe et al. 1999). Phytoplankton has a degradation index of 1 to 1.5 and a first-order rate 
constant of >10 yr’. Progressive degradation of organic matter lowers the degradation index and 
the reactivity constant 
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linked to organic matter composition. This does not imply that organic matter 
changes are the sole cause, but it does provide avenues to link kinetic biogeo- 
chemistry and organic geochemistry. 


Box 6.1: Ecological stoichiometry 
The traditional approach towards food-web and ecosystem functioning 
focuses on carbon, i.e. energy, flows and ignores the role of other elements 
and the biochemical makeup of organic matter. In other words, the traditional 
approach focuses only on quantity and largely ignores the quality of the 
organic matter. This is unfortunate, as all organisms need at least 20 elements 
or more and many organisms also have to acquire some essential compounds 
from the environment or their diet. Ecological stoichiometry explicitly deals 
with the flow of nutrients (primarily N, P, but also others) from the envi- 
ronments via primary producers to herbivores, detritivores and carnivores. 
Inspired by the pioneering work of Redfield in the early 20th century, 
ecologists (Sterner and Elser 2002) have developed the theory of ecological 
stoichiometry in which homeostatis is a central concept. Homeostasis is the 
property of organisms to maintain a constant composition, despite living in a 
variable environment or feeding on a variable diet (Fig. 6.7). Some organisms 
lack homeostasis, and their cellular nutrient levels reflect that in the envi- 
ronment (you are what you eat), while strict homeostatic consumers maintain 
their composition. Heterotrophs are, in general, more homeostatic than 
autotrophs, but intermediate behavior has been reported. Phytoplankton has 
been shown to have high flexibility in term nutrients, but suboptimal nutrient 
contents come at the expense of performance (e.g. growth rate). Consumers 
living on resources that deviate from their tissue composition have multiple 
ways of elemental adjustment: food selection, feeding behavior, regulation of 
assimilation and metabolism (e.g. diverting towards respiration or excretion). 
Most herbivores are richer in N and P than primary producers; this stoi- 
chiometric mismatch is smaller in aquatic systems than in macrophyte sys- 
tems in which carbon-rich lignin and structural cellulose are more abundant. 
Mismatches between the composition of a consumer and its resources has 
consequences for trophic transfer efficiencies and growth rates. The growth 
rate hypothesis involves a direct link between the growth rate, phosphorus 
content and ribosomes: i.e. high P content imply high RNA and thus high 
potential for growth. Ecological stoichiometry links the elements via com- 
pounds to the functioning of organisms and organism-scale processes to the 
global biogeochemical cycles. 
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Fig. 6.7 Generalised stoichiometric pattern relating consumer (N:P, C:N, etc.) ratio to resource 
ratios. Organisms lacking homeostasis reflect the ratio of their resource (dashed 1:1 line), while 
strictly homeostatic organisms maintain their composition irrespective of the environment or their 
diet (horizontal line) 
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